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SINGULARITY ANALYSIS FOR HEAVY-TAILED RANDOM 

VARIABLES 

NICHOLAS M. ERCOLANI, SABINE JANSEN, DANIEL UELTSCHI 


Abstract. We propose a novel complex-analytic method for sums of i.i.d. 
random variables that are heavy-tailed and integer-valued. The method com¬ 
bines singularity analysis, Lindelof integrals, and bivariate saddle points. As 
an application, we prove three theorems on precise large and moderate de¬ 
viations which provide a local variant of a result by S. V. Nagaev (1973). 
The theorems generalize five theorems by A. V. Nagaev (1968) on stretched 
exponential laws p{k) = cexp(—k^) and apply to logarithmic hazard func¬ 
tions cexp(—(logfc)^), jd > 2; they cover the big jump domain as well as the 
small steps domain. The analytic proof is complemented by clear probabilistic 
heuristics. Critical sequences are determined with a non-convex variational 
problem. 


1. Introduction 

The motivation of the present article is two-fold. First, we present a new analytic 
method for the investigation of large powers of generating functions of sequences 
that satisfy some analyticity and log-convexity conditions. The method is explained 
and developed for probability generating functions but it has potentially broader 
applications and is motivated by techniques commonly used in analytic combina¬ 
torics m- Specifically, we show that methods akin to singularity analysis can be 
pushed beyond the realm of functions amenable to singularity analysis in the sense 
of ffH Chapter VI. 1]. 

Second, we explore consequences for probabilistic limit laws and prove three the¬ 
orems on precise large and moderate deviations for sums of independent identically 
distributed (i.i.d.) random variables that are heavy-tailed [7] and integer-valued. 
The theorems generalize results on stretched exponential laws by A. V. Nagaev [TS] 
which have recently attracted interest in the context of the zero-range process [5]. 
They are close in spirit to results by S. V. Nagaev [m, however with more concrete 
conditions on the domain of validity of the theorems, and provide deviations results 
“on the whole axis” [20) . Our assumptions are more restrictive than one may wish 
from a probabilistic perspective; in retnrn, they allow for sharp results and may 
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provide a helpful class of explicit reference examples. For example, we prove that 
one of the bounds of the (local) big-jump domain for logarithmic hazard functions 
derived in [1] is sharp. 

The analytic proof of the theorems is complemented by clear probabilistic heuris¬ 
tics. Our results cover different regimes: a small steps or moderate deviations 
regime, where a classical variant of a local central limit theorem with corrections 
expressed with the Cramer series holds |12j . and a big-jump regime where the large 
deviation is realized by making one out of the n variables large. In the language 
of statistical mechanics and the zero-range process, they correspond to supersatu¬ 
rated gas and a condensed phase [2]. The critical scales that distinguish between 
regimes are defined with the help of a non-convex variational problem which en¬ 
codes competing probabilistic effects. This complements the approach taken by 
Denisov, Dieker and Shneer [3]. 

The study of combinatorial generating functions shares much in common with the 
study of probability generating functions; in fact in many instances they coincide 
or run parallel as is the case for more recent investigations in the area of random 
combinatorial structures [1]. From the viewpoint of complex function theory the 
key here involves relating asymptotic questions to questions about the nature of the 
singularities of generating functions viewed as more global analytic objects. In most 
of the successful applications of this singularity analysis to coefficient asymptotics 
a bridge is provided by the realization that the series in question satisfies some 
global algebraic or differential equation. Generating functions for which this is the 
case are referred to as holonomic. Pushing beyond this class in a systematic way 
requires new ideas and one of the most promising of these is the use of Lindelof 
integrals. Lindelof introduced these classically m as a means to constructively 
carry out analytic continuations of function elements (series) in a fairly general 
setting. In more recent times his construction has begun to be used to study 
non-holonomic combinatorial generating functions |10j . The generating functions 
for heavy-tailed distributions studied in this paper are of non-holonomic type and 
our methods of studying them show a new application of Lindelof’s construction 
that has novel connections to other areas of analytic asymptotic analysis such as 
bivariate steepest descent. In future work we hope to build on the present article 
in a way that broadens the application of harmonic analysis and complex function 
theory to problems of asymptotic analysis in both probability and combinatorics, 
such as applying the theory of Hardy spaces and Riemann-Hilbert analysis and 
extensions of Tauberian theorems as originally envisioned by Paley and Wiener [5] . 

Our proof shares some features with where cumulative distribution func¬ 
tions are approximate Laplace transforms and approximating moment generating 
function admit analytic extensions (see Section]^. Contour integrals that appear 
in inversion formulas are deformed and analyzed by Gaussian approximation—our 
proof details in Section 5.4 follow m- There are, however, key differences: we 
need not deal with approximation errors because of stronger analyticity assump¬ 
tions, and our detailed analysis of the underlying variational problem allows us to 
formulate more concrete conditions for our theorems. 

The remainder of this article is organized as follows. In Section we formu¬ 
late our main results and discuss applications to stretched exponential weights 
cexp(—fc“)) and weights cexp(—(loga:)^) with logarithmic hazard functions. In 
Section we explain the proof strategy in five steps, which are treated in detail 



SINGULARITY ANALYSIS FOR HEAVY-TAILED RANDOM VARIABLES 


3 


in the remaining sections. Steps 1 and 2 concern analytic extensions and notably 
use the Lindelof and Bromwich integrals (Section]^. Steps 4 and 5 analyze the 
critical points of a bivariate function and deal with the Gaussian approximation 
to a double integral (Section]^. The pivotal Step 3 connects the contour integral 
and the bivariate double integral; it leads to the full proof of our theorems that can 
be found in Section In Section we sketch how the method developed in this 
paper may be extended to more general settings, in particular when the analyticity 
required by our proof methods is not a priori given. 


2. Results 


We use the notation a„ ~ a„ = (1 + o(l))6n and a„ ^ bn = o{bn). 


2.1. Preliminaries. In order to formulate the results, we need to introduce crit¬ 
ical sequences deduced from a variational problem and the Cramer series. Let 
X, Xi, X 2 , ... be independent, identically distributed random variables with values 
in N and law 

P(X = k)= p{k) = exp(- 9 (fc)) {k G N) (2.1) 

for some sequence {q{k))ken- We assume that X is heavy-tailed and has moments 
of all orders, i.e., the generating function 

00 

G{z) = J2p{k)z'^ ( 1^1 < 1 ) ( 2 . 2 ) 

k=l 

has radius of convergence 1 and E[X’"] = k™‘p{k) < 00 for all m € N. Let 

p and cr^ be the expectation and variance of X. Set Sn = Xi + ■ ■ ■ + Xn- We 
are interested in the asymptotic behavior of P(S'„ = pn Nn) when n, Nn —> 00 
with Nn ^ y/n. The following assumption is similar to conditions considered by 
S. V. Nagaev [T7] . 


Assumption 2.1. For some a > 0, the sequence ( 9 (A:))fcgNn(o,oo) extends to a 
smooth function q : (a, 00 ) —>■ K which has the following properties: 

(i) q' > 0 , q" < 0 , and 9 ^^) > 0 . 

(ii) lima;_>oo a:g'(ai)/(loga;) = 00 . 

(hi) Ci^-^ < \q"{x)\ < C 2 ^-^ for some constants ci,C 2 > 0 . 

(iv) for some constants 03,04 > 0 . 

(v) q'{x) < for some a G ( 0 , 1 ). 


Assumption |2.1| allows for an easy analysis of an auxiliary variational problem, 
which is essential to the formulation of our main results. Let us collect a few 
elementary consequences. Under Assumption 2.1 q is concave on (a, 00 ) and p = 


exp(— 5 ) is log-convex. Moreover, lim, 
using 

q'iy) 


, x^q''(x)/ log X = —00 and ior y > x > a, 


we estimate 


q'{x) 

(I) 


= exp 


(-/* 


^ (fiyl ^ 

- q'{x) - 


q'{u) 


dit 


< 1 


(2.3) 


(2.4) 


Va;/ “ q"{x) ~ Vcc/ 


(2.5) 


Similarly, for y > x > a, 
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Since G{z) = J2k exp(—g(A:)) has radius of convergence 1, we also know that 


lim q'{x) = lim q"{x) = lim = 0 . 


( 2 . 6 ) 


Indeed by Assumption 2.1 q' is eventually decreasing and the limit i := lim^-oo q'{x) 
exists in M U {—oo}. Then i = lima;_>oo q{x)/x and G(z) has radius of convergence 
exp(£) = 1, whence ^ = 0. Assumptio n |2.1[ iii) and (iv) leads to the statements on 
higher order derivatives. Assumption |2.1[ v) implies q{x) = 0(a;“) as a: —> oo. 

Our method of proof requires two more analyticity assumptions, though in Sec¬ 
tion]^ we discuss how, following one might be able to relax these assumptions. 

Assumption 2.2. There exists & > 0 such that (p(n))„gpjn[b,oo) extends to a func¬ 
tion p{^) that is continuous on a closed half-plane Re^ > b, analytic on the open 
half-plane Re^ > 5, and in addition satisfies 

(i) For every e S (0,7r), some G^ > 0, and all we have |p(^)| < exp(e|^|). 

(ii) 1(6 -I- \s)^p{h + is)|ds < oo for all A: € N. 


Moreover p(a;) = exp(—( 7 (x)) for all x > max(a, 6 ) with a, q{x) as in Assumption 2.1 


Assumption 2.3. Let p{^) = exp(—g(^)) be the analytic extension from Assump¬ 
tion 2.2 defined in Re^ > 6 . Then q{^) = —Log^, defined with the principal branch 


of the logarithm is analytic as well, and the following holds: 
(i) For r > 0 large, let = b + — 6 ^. Then as r —>■ 


'ReC=b. |5|>r 


exp(-Reg(0)d^ 


< exp(—Reg(zj.) + O(logr)). 


(ii) Im (^q'(^)) < Im {^q'{r)) for all large r and all ^ with Im^ > 0 and | 5 | = r. 

(iii) \q^^\0\ < C'|g"(C)/CI for some C > 0 and all 

Assumption |2.3| enters the proof of Theorem |4.4| only. 


Variational problem and critical scale. Assumption |2.1| is tailored to the analysis 
of an auxiliary variational problem, motivated by the following heuristics. For 
subexponential random variables, the typical large deviations behavior is realized 
by making one out of the n variables large, 

V{Sn = tJ.n + Nn) « n¥{Xn = Nn- kn)V{Sn-l = fin + kn) (2.7) 

with a yet to be determined optimal k„. Assuming that a normal approximation 
for the second factor is justified, we get 

F{Sn = pn + Nn) « exp(^-q{Nn - kn) - (2-8) 

where we have neglected prefactors n and (see Eq. (2.11) below for a 


more refined heuristics). The optimal kn is then determined by minimizing the 
term in the exponential. Thus we are led to the minimization of 

{Nn -X)^ 


fnix) = q{x) + 


2na‘^ 


(2.9) 


As illustrated in Figure the nature of the variational problem changes with A„. 
Define x* > 0 and N* by 

1 


9"«) = - 


n: = '. 


ncr^ 




( 2 . 10 ) 
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For sufficiently large n, the inflection point x* is uniquely defined because of the 
monotonicity from Assumption 2.1 and Eq. (2.6), moreover cc* —)■ oo. The quantity 
N* is defined in such a way that the tangent to the curve y = q'{x) at x = x* has 
equation y = {N* — x)/(n(T^), see Fig. 

The next two lemmas characterize the minimization of /„; they are proven in 
Section [O] The Hrst lemma relates the critical points of /„ to the location of Nn 
compared to N*. 


Lemma 2.4. For sufficiently large n, the following holds true: 

(a) If Nn < N*, then /^ > 0 on (a,oo). 

(h) If Nn = N*, then fn has the unique zero x*, moreover ffix) > 0 with 
equality if and only if x = Xn- 

(c) If Nn > N* and limswpn^oo < lim 2 ,\,o then ff has exactly 

two critical points x„ and x'n, which satisfy x'n < Xn < Xn < Nn and 

fn {x'n) = max fn , fn{Xn) = ^ min fn . 

(a,tc*) (x*,oo) 

(d) If Nn > N* ond lim inf„_).oo iV„/(n(T^) > lima;^^^ q'{x), then fn has a unique 
critical point x„. It satisfies x„ € (x*, A^„) and is a global minimizer. 

For Nn > N*, the function /„ may have two local minimizers: a and x„, and we 
may wonder which one is the global minimizer. The answer depends on the location 
of Nn compared to a new critical sequence N**- Concrete examples are given in 
Sections 12..31 and 12.41 


Lemma 2.5. For n sufficiently large, there is a uniquely defined N** > N* such 
that: 

(a) If N* < Nn < N**, then /„(a) < /„(x„). 

(b) If Nn = N**, then /„(a) = /„(x„). 

(c) If Nn > N**, then /„(x„) < /„(a). 

In general it may not be straightforward to determine N** exactly, but it is simple 
to find a lower bound: if a = 0 and q{Nn) < A^^/(2ncr^), then Nn > N**. This 
lower bound corresponds, roughly, to the sequence A(n) in [T7] . 

The sequences introduced up to now are ordered as follows. 

Lemma 2.6. As n —>■ oo, we have 

ffn <s: x*n < N* < N** = ^ n, 

and for some constants C,S > 0 

{1 + 6)x*n < K < Cx*n. 


The lemma is proven in Section [O] Lemmas 2A and |2.5[ together with the heuris¬ 
tics described above, suggest that for Nn > N** the unlikely event Sn = ny + Nn 
is realized by making one component of the order of x„. One may wonder how far 
Xn is from swallowing all of the overshoot Nn- 

Lemma 2.7. Suppose Assumption holds true. Let Nn > N* ■ Then 

' N* 


Nn-K<Xn^ 


< Nn, na^fn{xn) = 1 “ |(?"(x„)| = 1 + 
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Figure 1. Minimization of fnix) = q{x) + {Nn — a;)^/(2ncr^) and 
illustration of Lemmas 2.4 and|2.5|for weights q : (0, oo) —)■ M with 


g(0) = 0. For 7V„ > N^, fn{x) has two critical points x'^ and 
Xn separated by an inflection point a;*. The global minimum is 
reached either at a: = or at a; = 0. 


In particular, for iV„ ^ iV*, we have ^ Nn and na^f''{xn) —>■ 1- The lemma 
is proven in Section |5.1[ The information on the second derivative enters a refined 
heuristics: we make the ansatz that conditional on the unlikely event S'„ = nfj,+Nn, 
there is one large component of size a;„, but the size is not deterministic. Instead, 
there are fluctuations around This yields 


P(S'„ = n/a + Nn) ~ ^ = Xn+ £)P(5'„_1 = Nn - Xn - ■ 

e 

exp(-/„(a:„ +£)) 




V2 


irncr^ 


n e 


-fniXri) 


E 


exp(-/"(a:„)fV2) 
V 27rncr^ 


_ nexp(-/„(a:„)) 

a /1 - na'^\q''{Xn)\ 

Theorem |2.11| below confirms the heuristics for large Nn, up to correction terms 
both in the prefactor and in the exponential. 


Cumulants and Cramer series. The heuristics together with Lemma 2.7 suggest 
that the optimal kn = Nn — in Eq. (2.8) is of order up to N* ^ a/u. At this 


scale the normal approximation fails and requires correction terms. The latter are 
usually expressed with the Cramer series m, whose definition we briefly recall. 
Let ip(t) be the cumulant generating function of X, 


= logE[e‘^] =logG(e*) (Ret < 0). 

Note :/3(0) =0. As t —>■ 0, (f{t) can be approximated to arbitrary order 


( 2 . 12 ) 


(2.13) 

i=i 

with finite and real expansion coefficients kj G K, the cumulants, see Section 
Notice that ki = /i is the expectation and K 2 = > 0 is the variance of X. 
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Definition 2.8. Let tir) = -\- J2j>2 UjT^ be the formal power series obtained 

by inverting 


J>2 


j'- 


The CramCT series J2j>o ^^6 formal power series defined by composing the 

expansion of t(T) with the expansion of (/i -|- r)t — tp{t), 




i>i 


t{T} 


j>0 


Equivalently, the Cramer series is the left-sided Taylor expansion of the Legendre 
transform (p* at fi: let p*(x) := sup(<;g(ta: — <p{t)). Then as r ^ 0, 


P*{p + t) = - 


2 a 2 


-I- ^ XjT^ + 0 {t 

7=0 


r+4 


(2.14) 


to arbitrarily high order r. 


Remark 1. For t > 0, log[X)fc>i Pik) exp(fct)] is infinite and the standard convention 
is to set (p{t) = oo; then (p*{yL+T) = 0 for t > 0 and Eq. (2.14) no longer applies. We 
adopt a different convention, however, for which ip{t) is smooth in a neighborhood 
of 0 (see Theorem 4.2), though it becomes complex-valued, and Eq. (2.14) applies 
to (Re ip)* for positive r as well. 


2.2. Main theorems. Set fno = fn and for r > 1, 

(7V„ - xf (7V„ - a;)^ 


fnr{x) = q{x) -I- 


2 ncr 2 


r—i 


Nr,. - X\i 


(2.15) 


7=0 


Remember the minimization of fnM and the critic al sca les i/n <C iV* < N** = 
0 (n^/’^ 2 -a)) ^ illustrated in Fig. Ill In Proposition 


2.12 


below we check that the 
properties of f„ carry over to fnr- The following theorems provide a local variant 
of a large deviations theorem by S. V. Nagaev see also [m Theorem 2.1]. 
The principal difference, apart from the local character of the theorems, is that our 
detailed investigation of the variational problem and notably Lemma |2.7| allows us 
to formulate conditions directly in terms of Nn, whereas S. V. Nagaev’s criteria 
included an indirect condition on the sign of some second derivative. 

Theorem 2.9. Let Nn — )■ oo with ^/n <C Nn < (1 + o{l))N*. Pick r large enough 
so that n{Nn/nY —>■ 0. Then 

iV 2 

*71 _|_ 

-.2 


P(5'„ = pnP Nn) 


1 




:exp 


2 ncr 2 




7=0 


Theorem 2.10. Let Nn —t oo with liminf Nn/N* > 1 and Nn = “1). Pick 

r large enough so that n(7V*/n)’’ —> 0. Then 


P(S'„ = pn + Nn) = (1 + o(l)) 


1 


: exp 


-(1 + 0 ( 1 )) 


/ N^ Ny^ cNny\ 

V 277(72 + ^ ) ) 

^ fnr {Xnr)'^ 


277(t2 
77 


7=0 


- na'^\q"{Xnr)\ 


exp -J 
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with Xnr = Nn + 0{N*) the largest solution of f^^{xnr) = 0. 


Lemma 2.5 suggests that for Nn ^ ^^e first contribution dominates and for 

Nn ^ Nf* the second contribution wins, but one has to be careful because of the 
factors n and l/\/27rn(T^ as well as the Cramer corrections; a detailed evaluation is 
best left to concrete examples (see, however. Corollary 2.13 below). 


Theorem 2.11. Let Nn 

n{N*/nY 0. Then 


00 with Nn S> . Pick r large enough so that 


P(5'„ = iJ,n + Nn) ~ nexp^-/„r(a:nr)) 


with Xnr = Nn + 0{Nf) the largest solution of fn^i^nr) = 0 . 

In practice one may prefer not to deal with the Cramer corrections or the variational 
problem, and the following proposition is helpful. 

Proposition 2.12. Suppose that liminf„_,,oo 1V„/7V* > 1. Fix r S Ng. Then, 
for sufficiently large n, Xnr is the maximizer of fnr restricted to (a;*, and the 
unique zero of f^r in that interval. Moreover 1 — n(T^|q"(x„r)| = 1 + 0{Nn/Nn) 
stays bounded away from zero and 

Xnr = iv„ - (1 + o{l))na^q'{xnr) = Nn + O^N*), 
fnr(Xnr) = q{Xnr) + ^(^ + o{l))na‘^ q'{x nr) = + ‘^(^)) ' 

The proposition is proven in Section 15.11 For iV„ N*, we obtain Xnr Nn, 


(2.16) 


qfxnr) ^ qfNn), and q{xnr) = q{Nn) - (1 + o{l))na'^q'{Nn), hence 

fnr{Xnr) = q{.Nn) “ ^(1 + o{\))na‘^ q {Nn)'^ ■ 

Now suppose in addition that liminf 2 rnT^ lq{i^n) > 1- Then in Theorem 


2.10 


the 


first summand is of order exp(—(1+ o(l))iV,^/(2nCT^), the second of order exp(—(1 


o(l))g(iV„)), so the first contribution is negligible and the validity of Theorem 2.11 
extends accordingly, since 1 — na'^\q"{xnr)\ 1 for Nn N*. Eq. (2.161 now 
yields the following corollary. 


Corollary 2.13. Take Nn —>■ oo with Nn ^ N* and liminf /q{Nn) > 1- Then 
P(S'„ = npi + Nn) = n¥{X = Nn) 

if and only if 


In concrete examples. Theorem 2.9 should allow us to extend the domain of validity 


of the corollary to 1V„ ^ Nf*. The condition 


0 is closely related to 


the insensitivity scale discussed by Denisov, Dicker and Shneer [4], as 

p{Nn ± vW2) 


P{Nn) 


—> 1 y na'^q'{Nn) —>• 0. 


(2.17) 


Remark 2 (Big-jump vs small steps). The domain where P(S'„ = n/i -I- Nn) ~ 
nF{X = Nn) is sometimes called big-jump domain. Think of Sn as the position of a 


random walker with step size distribution p{k). In the situation of Corollary 2.13 
the unlikely event that the walker has travelled a distance /rn -I- Nn much larger 
than the expected distance /in is realized by one big jump of size iV„. Finding the 
boundary of the big-jump domain is an active field of research [3]. 

















SINGULARITY ANALYSIS FOR HEAVY-TAILED RANDOM VARIABLES 


9 


The interpretation of Theorem |2.9[ in contrast, is that the moderate overshoot 
iV„ is achieved by a collective effort: all steps tend to stay small, though each 
stretches a little beyond its expected value /r. For stretched exponential variables, 
this interpretation is made rigorous in |15j and [5]. 


Remark 3 (Condensation in the zero-range process). In the zero-range process, 
the random variables Xi,..., model the number of particles at lattice sites 
j = with Sn the total number of particles, and ^ is a critical density. 

Theorem |2.9| corresponds to supersaturated gas. In Corollary |2.13[ the particle 
excess Nn is absorbed by a condensate, i.e., one large occupation number. In 
Theorem 2.II the particle excess is shared by a condensate of size Xnr < Nn and 
supersatured gas. See [2] and [8l Section 7]. 

We conclude with an equivalent but more intrinsic formulation of Theorem |2.11| 
In Section we shall see that G{z) extends to a function that is analytic in the 
slit plane C \ [I,oo), and in addition the limit G(e‘) = lime\^o G(e‘ -I- ie) exists 
for all t > 0. So the cumulant generating function ip{t) = logG(e‘) extends to a 
function that is well-defined and smooth in a neighborhood of the origin in K, and 
Eq. (2.131 stays valid for small positive t. We define 

^n{t, C) = -q{0 + fit) - + Nn)t + 


•p _ 

= -q{0 -t{Nn-0 + 0{f). 

3=2 


(2.18) 


The asymptotic expansion holds for every order r. For Nn ^ and sufficiently 
large n, the bivariate function $„(t, 5 ) has exactly two critical points in ( 0 , ^^) x 
(a, oo). We label them as (t„,C„) and {t'n,^'n) with < t'n- Then, in the situation 
of Theorem |2.1H we have 

P(S'„ = n/r -f iV„) ~ n . (2.19) 


It is in this form that we prove the theorem. Let us explain how to recover the 
expression in terms of fnr- We may solve for V<i>„(<,^) = 0 in two steps: (1) use 
dt^n{t,C) = 0 to express t = t{^) as a function of ( 2 ) plug the expression into 
9^<I>„(t,^) to obtain an equation for This latter step breaks into the following 
two stages: 

( 2 a) substitute the expression of t = t{^) into the expression for $„(t,^) so as 
to obtain a function ‘&n(t(C ))0 of C only; 

( 2 b) set the derivative of with respect to ^ to zero. 


which is valid since 

(i fif 3 ^ 

By the definition of the Cramer series, step (2a) gives 

{Nn-0^ , (Nn-O^^, /Nn-^y 


= -Qio - 


2ny 




3>0 


( 2 . 20 ) 


( 2 . 21 ) 


Truncation of the asymptotic expansion on the right-hand side gives precisely the 
function —fnriO- For r large enough, step ( 2 b) shows that we may approximate 

^nitn,^n) = -fnriXnr) + o(l), na^q''{^n) = na'^q'\Xnr) + o(l) ( 2 . 22 ) 


hence the equivalence of Eq. (2.191 with the expression from Theorem 2.11 
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2.3. Application to stretched exponential laws. Here we explain how to re¬ 
cover five theorems by A. V. Nagaev [TS] for stretched exponential variables. Let 
a G (0,1), c > 0, and 


p{k) = cexp{—k°‘), q(fc) = fc“ —logc (fc G N). (2.23) 

We need not check Assumption |2.3| since Theorem |4.4| for stretched exponential 
weights has already been proven in [121 Theorem 2.4.6]. 


Lemma 2.14. The probability weights (2.23) satisfy Assumptions 2.1 and 2.2 and 
we have 


= [q;(1 — a)na‘^] 


2ll/{2-a) 


1 — a 


with Ca = (2 — a)(2 — 2a) “)/(2 Moreover Vn^qfNn) 0 if and only if 

Nr, » 


The proof of the lemma is sketched in Appendix [a] The critical scale n is 

explained by a simple scaling relation; for we have 






(2.24) 


A careful examination of the expressions in Theorem |2.10| shows that the first 
summand dominates if < (1 — 6)Nf^* while the second dominates if A„ > 
(1 -I- (5) A** for some 5 > 0. In particular, Theorem 12. 9| e xtends to A„ < (1 — S)N**, 
which corresponds to Theorem 1 in [15]. Theorem 2.10 for A„ ^ A** is Theorem 4 
in [T5|. For A„ > (1 -|- <5)A**, we have 

1 


P(S'„ = np + Nn) 


3-Ar(a:,ir-) 


y'l - na^q'fxr 


0 


(2.25) 


This regime can be divided into three cases: 

(a) Nn corresponds to Theorem 2 in [15j . 

(b) When A„ is of the order of the corrections from the Cramer series 

are irrelevant and 


P(S'„ = pn + Nn) ~ nexp(-/„(x„)), (2.26) 

i.e., we may choose r = 0. This corresponds to Theorem 6 from the erratum [IBl, 
replacing Theorem 5 in the original article m- The statement actually extends 
to Nn S> Indeed q'{xnr) (Nn — Xnr)f{na'^) = o{Nn/n) yields 

Nn — Xnr = (1 + o(l))ncr^aA““^ and 

fnriXnr) = fn{Xnr) + O(uA") , (2.27) 


2.13 


and one can check that fnixnr) = fn{Xn) + o(l) if A„ n 3“). 

(c) Nn ^ 7 T,“i/( 2 - 2 a) corresponds to Theorem 3 in [T5| and our Corollary 

2.4. Application to logarithmic hazard functions. Here we specialize to 

p(fc) = cexp(—(log/c)^), 9 (fc) = — logc-b (log/c)^ (2.28) 

with P > 2[^ 


^For /3 e (1, 2 ], xq'{x) = ,S(log3;)^ oo but the stronger condition xq'{x)l logic oo from 

Assumption |2.l| ii) fails. We suspect that this restriction is technical and could be lifted with 
more detailed estimates, but a proof or disproof is beyond this article’s scope. 
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Lemma 2.15. The weights ( |2.28 ) satisfy Assumptions 2.1- 2.3 Moreover 
N* - 2x1 - 2^2i-/5/3riCT2(iogn)/3-i, Nf* ~ 2na'^[\ognY, 

and 'Jncr^q'{Nn) —)■ 0 */ and only if Nn ^ y/n{\ogn)^~^. 

The lemma is proven Appendix[S] Notice that unlike the stretched exponential case 
(Lemma |2T4t , Nff is much larger than . The scaling relation ( |2.24[ ) is modified 
as follows: for iV„ = k^Jna‘^{\ogn)P~^ ^ kx'^, we have 


fniyxl) = {logx*jf^ + Pilognf ^(\ogy 


(k-yf 


) +o((logn)'^-i) (2.29) 


and (logx*)^ -- (logn)^ -- (logiV„)^. 

Our results may now be applied to obtain a sharp boundary for the big-jump 
domain. 


Theorem 2.16. Let p{k) be as in Eq. (2.28) with (3 > 2 and Nn ^ y/n. Then 
]P(<S'n = yn + Nn) ~ n¥{X = Nn) if and only if Nn ^/n{logn)^~^. 

The “if” part of the theorem is actually a special case of [U Theorem 8.2] and as 
such not new. The “only if” part shows that the boundary derived in [3] is in fact 
sharp. 

Proof Theorem 2.16 Let /„ := yJna'^{\ogn)^~^ and notice In ^ N** for /3 > 2. 
Suppose that iV„ J„. Then we have, in particular, Nn Nf^*. Write Nn = 
anNff with a„ —)■ oo. Then 


(logiV„)/3 




logn 


> 


(loga„)/5 


(2.30) 


and it follows from Corollary 2.13| that P(5'„ = ny + Nn) ~ nF{X = n), hence 
Nn ^ N** is indeed a sufficient condition. In order to check that it is necessary. 


we treat the case = 0(/„) with Theorems 2.9 and 2.10 


Case 1: N* Nn = 0{ln). In Theorem 2.10 we obtain a lower bound by 
neglecting the first contribution and estimating 1 — na‘^\q"{xnr)\ < 1- Combining 
with Proposition |2.12[ we find 

P(S'„ = nyL + Nn) > neJip{-fnr{Xnr) + o(l)) 

> nex.p(^-q{Nn) + (1 + o(l))n(T^g'(fV„)^ -f o(l)^ (2.31) 

hence in view of Nn = 0[ln) and Lemma[2.15| 


lim inf 


P(S'„ = n/r -b Nn) 


> liminf exp((l-b o(l)ncr^g'(A^„)^)) >1. (2.32) 


nF{X = Nn) 

Case 2: Nn = 0{Nff). Then we have in particular Nn = o{N**). Write iV„ = 


UnNff with 


0, then 

K < a^(logn)^ 
2ncr2(logiV„)/3 “ (logv^n)^ 


and Theorems 12.91 and 12.101 show 

P(5'„ = ny + Nn) Nf 
log -- > - 


nP(X = Nn) 


2710-2 


(log Nn)^ + 0(log n) ^ OO. 


(2.33) 

(2.34) 

□ 
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3. Proof strategy 


Here we explain the strategy for the proof of Theorems |2.9f|2.1l] We focus on 
the case Nn = o(n) and Theorem 2.10 Set m = fj,n + Nn- We start from the 
observation that P(S'„ = m) is equal to [z™] G'(z)™, the coefficient of z™ in the 
expansion of G{z)'^, which in turn is given by contour integrals 


/•27ri 

I dt. (3.1) 

0 

The contour integral can be taken over any circle centered at the origin with 
radius r < 1. A steepest descent ansatz would look for a point Zn such that 
ZnG'{zn) = m/n = p. + Nn/n, or with ip'(rin) = /r + Nn/n, and then integrate 
over \z\ = \zn\ (or Ret = Rery^). However in the regime m/n > G"(l) = yi that 
we investigate there is no such point, and instead we follow an approach that is in 
the spirit of singularity analysis m but with several novel ingredients. Crucially, 
the generating function G{z) does not fall into the class of functions which Flajolet 
and Sedgwick call “amenable to singularity analysis” [TT] [Chapter VI]. 


[z™]G(z)" = 


27ri 


G(z)” dz _ 1 
z™ z 27ri 


Step 1: Analytic extension to slit plane. Observe that G(z) has an analytic 
extension to the slit plane C\[l,oo). This is proven with the help of the Lindelof 
integral [nmni, see Proposition |4.1[ The key ingredient here is that p(^) is ana¬ 
lytic in a complex half-plane containing the integers fc € N and growth slower than 
exp((7r — e)|^|) as ^ oo. 


Step 2: Behavior near the dominant singularity and along the slit. 

The Lindelof integral actually shows that the analytic extension G(z) has well- 
defined limits as z approaches the slit [l,oo) from above or below, i.e., the limits 
lim^^o G{ e* -l-ie) and lim^x^o e* — ie) exist for all t S K. Moreover the imaginary 
part along the slit is given by a Bromwich integral, 

1 ^l/2-l-ioo 

lim ImG(e* -I-ie) = — / e*^p(^)d^ (1 > 0). (3.2) 

2l Ji/2_ioo 

The line of integration Re^ = 1/2 can be replaced by any other line Re^ = x > 0. 
As t \ 0, the imaginary part vanishes faster than any power of t, whereas the real 
part can be approximated to arbitrarily high order by a Taylor polynomial. 


Step 3: Contour integrals. We may now deform the contour of integration: 
in the z-plane, we replace the circle of radius 1 by a Hankel-type contour consisting 
of a circle of radius e’^ and a piece hugging the segment [l,e), see FigureIn the 
t-plane, we replace the vertical segment joining 0 and 27ri by the three other sides 
of the rectangle with corners 0, e, e -I- 27ri, 27ri. This yields 


1 


r2Tr 


[z’"]G(z)” = / e^v{e+i9)-m{e+W) 

27ri \Jq Jq 


^uv{t+2ni)-m{t+27ri) ^3 3 ) 


We focus on = o{n) and choose e = ijn as the solution of 


Rev? (»7n) = — = fJ.-\ - 

n n 


(3.4) 
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(b) 


Figure 2. Contour integrals in the z-plane (a) and in the t-plane 
(b). Dotted lines: Eq. (3.1). Solid lines: deformed contours in 
Eq. (3.3). Recall that z = e*. Later e = ??„ will be chosen in a 
judicious way. For Nn <C fV**, the dominant contribution should 
come from the horizontal pieces of the deformed contour in the 
t-plane. For ^ fV**, the dominant contribution should instead 
be from the vertical line. 


Notice r]n ~ Nn/{na'^). Using the identities 

G{z)=G{z), Vt > 0 : Lp{t + 2T:i) = ip{t), 


Eq. (3.3) becomes 


with 


= 


[z™]G(z)" = + U„ 


gnRe<^(t)-mt sin{nlm ip{t)) dt 


(3.5) 

(3.6) 

(3.7) 


Vn=^ [ cos(nImv5(r;„ + i0))d0. 

Jo 

Standard arguments show that the dominant contribution to Vn come from small 
9. Since ImG( e*) —>■ 0 faster than any power of t as t \ 0 and 

Im G( e*) = Im = e^® Im (p{t) ~ Im V5(t), (3.8) 

we may drop the trigonometric functions from Eq. (13.71) and find 


Hn-- 

TT 

TT 


e"R"‘^(‘)-'"‘lmG(e‘)dt, 

^nRc ip(r]„+ie)-mri„ 


(3.9) 


The vertical contribution is evaluated with the help of a Gaussian approximation 
around 0 = 0, which yields 


Vr,. 


1 


^J2^:nRelp"{r|n) 


j^nRe ^p{rin)—m7]ri 


(3.10) 
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With Definition |2.8[ we recognize in Eq. ( |3.10[ ) the asymptotic expression from 


Theorem 12.91 and obtain 


1 




: exp 


nna^ 




r—1 


2ntT^ 


j=o 


n 




(3.11) 


The evaluation of Hn is more involved. As a preliminary step, we express lmG'(e*) 
through the Bromwich integral and find 

]dt (3.12) 


Hr, 


n 

27ri 


r'lr, / ri/2+ioo 

1 Jo ^Jl/ 2 -ioo ^ 


with il>„ as in Eq. (2.18). 


Step 4: Critical points of $„(t,^). In order to apply a Gaussian approxima¬ 
tion to the bivariate integral ( |3.12 1, we look for a critical points (t„,^n) of with 
tn € (0,77„) and S (0, oo). The gradient V<i)„(t„,^„) vanishes if and only if 

tn — Q i^n), 

(3.13) 


Re - V5'(0)) = ——— 


Since Re(/?'(t) = fJ, + cr^t + O(t^) as f \ 0, Eq. ( |3.13 ) implies 

N - e 

//> \ sn 

^ - 


(3.14) 


Hess4>„(t,C) = 


(3.15) 


(3.16) 


We recognize the equation for the critical points of /„. Lemma [2.4| suggests the 
following: for Nn ^ N*, there should be no critical point, for N* Nn n, 
there should be two. Let us focus on the latter case and label the critical points 
as (<n,^n) and with In view of Lemmas 2.4 and 2.7 we expect 

« Xn and ~ 2 :(j, hence 

tn - q'{Nn), ^n = Nn+ 0{N*) 
and ^'n < Xn < ^n- The Hessian of <!)„ is 

^nReip''{t) 1 
1 

Using again Lemma [2.7| and ~ Xn, we expect 

detHess$„(t„,^„) =-1 - (l-f o(l))ncr^g"(^„) =-1-b <0 (3.17) 

thus is a saddle point. (More precisely, it is a saddle point of Re$„, but 

the abuse of terminology is natural and not problematic in our context.) On the 
other hand ^n ^ x'^ < x^ with 1 + na^q'{Xn) = 0 by definition of a:*, so we expect 

detHess$„(t(j,^(j) = -1 - {l + o{l))na^q''> 0. (3.18) 

Step 5: Gaussian approximation for In order to evaluate the double 
integral in Eq. (3.12), we use a good change of variables and a Gaussian approx¬ 
imation. Let ^{t) be the solution of q'{^) = t, so that 9^$„(t,^) = 0 if and only 
if ^ = ^{t). It is convenient to deform the contour and integrate along Re^ = ^{t) 
instead of Re^ = 1/2. The integral becomes 

pr)ri , pOO 

Hn-^J (y e'^"(‘’«(‘)+‘*)ds)dL (3.19) 
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A straightforward computation shows that Fn{t, s) = Sn{t,^{t) + is), a function of 
two real variables t, s, has a critical point at (t„, 0) with positive definite Hessian 


HessF„(t„,0) - “ 


det Hess$n(tn,^n) 


(3.20) 


see Lemma B.l F„ : (0,ry„) x K —>■ C has another critical point at (t(j,0), with 
negative determinant of the Hessian; later we show that it does not contribute to the 
integral. The evaluation of is concluded by replacing the double integral (3.19) 
by the integral of the Gaussian approximation around (t„, 0) which yields 




( 271)2 




27r y det Hess F„(t„,0) 


a/I det Hess$„(t„,^„)| 




The argument leading leading to Eq. (2.22) and Eq. (3.17) show 
„ -^ 


y^l - na^q"{xnr) 


(3.21) 

(3.22) 


4. Analytic continuation. Lindelof and Bromwich integrals 


Here we take care of steps 1 and 2, starting from Assumption |2.2[ For concrete¬ 
ness’ sake we write down the results for b = 1/2; they apply for general b with 
straight-forward modifications. Define 

1 /■1/2-l-ioo 

H'w) = - 7 ^ piOw^- -(w G C\(-oo,0]), (4.1) 

xTi 7i/2-ioo sm7r4 

the Lindelof integral with symbol p(^). 

Proposition 4.1 (|13)1. 

(a) A{w) is analytic in the slit plane C\(—oo,—1]. 

(b) In the unit disk |w| < 1, 

OO 

A(u;) = ^p(fc)(—re)*’ = G{—w). 

fc=i 


A detailed proof and many additional properties of A{w) can be found in |10j . 
Proposition 4.1 shows right away that G{z) = A(—z) has an analytic continuation 
from the unit disk to the open slit plane C\[l, oo). We use the same letter G{z) for 
the analytic continuation, and set (p{t) = Log G{ e*) with Log the principal branch 
of the logarithm, and Imt G [0, 27r). We prove the following additional properties 
of G{z) and (p{t). 


Theorem 4.2. 

(a) The boundary value G( e‘) = limex^o G( e* -I- ie) exists for all t gM. and is 
a smooth function o/t G K. 


(b) The imaginary part ImG( e*), t >Q is given by the Bromwich integral (3.21. 

(c) ipit) is well-defined and smooth in a neighbourhood of the origin; the deriva¬ 
tives Kj = are real. As t —)■ 0 in the strip Imt G [0, 27r), we have 

p{t) = LogG(e‘) = ^Kj-^+0{G+^). 

1=1 


to arbitrarily high order r. 
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In (c) z = e* is allowed to approach the slit [l,oo) as fast as we like; we may 
even take t real. Because the coefficients Kj are real, we find in particular that 
ImG'( e*) vanishes faster than any power of t as t —>■ 0, t S M. 


Proof of Theorem 4-2 For u G C in the closed strip Imu S [—7r,7r], define 




r*l/2+ioo 






'i/2-ioo sinTrC 


de 


(4.2) 


When Imu is in the open strip Imu € (—tt, tt), we have lu = e“ S C\(—oo,0] and 
T(u) = A(e“ ). Along the vertical line Re^ = 1/2, we have 


exp(^u) 


sin(7r^) 


exp(-«ImO r../, = i + 


exp(7rs) + exp(—TTs) 


(4.3) 


By Assumption 2.2'ii), since ^V(C) is integrable along Re^ = 1/2. Eq. (4.3) 


then shows that the integral defining L(u) is absolutely convergent, and it stays 
absolutely convergent if we replace the symbol p(^) by C^p(C)- Standard arguments 
for parameter-dependent integrals then show that L(u) is continuous on the closed 
strip, differentiable in the open strip, and we may exchange differentiation, limits, 
and integration, which shows that the restriction of L to the boundaries Im u = zLtt 
yield smooth functions. 

When z —)■ e* S [1, oo) along Im z > 0, we have lu = —z —>■ — e‘ along Im w < 0. 
Thus we may write w = e“ with Reu —>■ t and Imu = argic —tt. Therefore 


/ OO 

p{\+ is)' 

-oo 


.(l/2+is)t _ 


exp(7rs) 


exp(7rs) -I- exp(—TTs) 


ds. 

(4.4) 

This proves the existence of the limit and, in view of the above mentioned properties 
of L(u), the smoothness as a function of t. The complex conjugate is 


/ OO 

p(i _ is) 

-OO 


exp(7rs) 


exp(7rs) -I- exp(—TTs) 

/ OO 

p(i+is)eF/2+-F 

-OO 


ds 


exp( —TTs) 


exp(7rs) -I- exp(—TTs) 


ds. (4.5) 


Therefore 


-I poo -I ^1/2+ioo 

ImG(e*) = -/ p(i+is)e(i/2+-)*ds=- / p(e)e*«de 

J— OO «/l/2 —ioo 


This proves (b). For (c), consider first real t G R. We have already checked (a) 
hence G( e*) is in G°°(M). It is real and strictly positive for t < 0 (this follows from 
the series representation and p{k) > 0), and non-zero though possibly complex¬ 
valued for sufficiently small t > 0. Therefore (p{t) = log G( e*) is well-defined and 
smooth in some interval (—oo,5), 5 > 0, and real-valued for t < 0. In particular, 
the derivatives kj = (0) exist and are real, and pit) can be approximated to 

arbitrarily high order by Taylor polynomials. The extension to complex t, Imt G 
[0,27r), follows again from the smoothness of T(u) in the closed strip Imu G [— tt, tt]. 

□ 
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Theorem |4.2[ b) has an interesting consequence. Eq. 
the formula for the inverse Laplace transform, therefore 

poo 

p{X) = TT / e“*'’'Im G( e‘)dt (ReA>0). (4.6) 

Jo 

In the special case of stretched exponential weights, we can draw on an extensive 
literature as exp(—A“) is known to be the Laplace transform of a probability density, 
an a-stable law. For a = 1/2 [S] 


(3.2) is, up to a factor tt. 


ImG(e‘) = |0e-^/(4*) (t > 0). 

For general a G (0,1), we have instead [T^l Theorem 2.4.6] 

~exp(-(l - 


ImG(e') 



t(2-a)/(2-2a) 


(4.7) 


(4.8) 


2.1 


as < \ 0. This is proven in [12] by applying a steepest descent approach to the 
Bromwich integral. For general weights, Eq. ( |4.8| ) is generalized as follows. 

Assume that t < lim 2 :\,a W'{x)\. By Assumption 
and negative on (a, oo). By Eq. (2.6), we have q''{x) - 
there exists a uniquely defined ^{t) that solves q'{^{t)) = t. We define 

= t^{t) - q{^it)) 

and note the relations 

1 


q” is strictly increasing 
0 as a; —)■ 00 . Consequently 


(4.9) 


= rit) = 


q"{my 


(4.10) 


so Tpy) is monotone increasing and strictly concave. Since is the Legendre 

transform of the convex function —q(x), it comes as no surprise that Assumption |2.1| 
on large x translates into information on small t. 

Lemma 4.3. The following holds: 

(a) limtx^o (i) / log t = 00. 

(b) for some c > 0 and all sufficiently small t > 0. 

(c) 0 < < g 1 ^^ for some G > 0 and all sufficiently small t > 0. 

The lemma has been proven in m Lemma 2.2]. 


Theorem 4.4. As t \ 0, 

ImG(e‘) ~ ^y/27rlV'"(t)| . 

Proof. By Theorem |4.2[ b) we may start from the Bromwich representation of 
ImG(e‘). The analyticity of q{£,) allows us to replace the contour Re^ = 1/2 
byT = riUr 2 where 

Ti = {e G C 1 Id = m, Ree > 1/2}, r2 = U G C 1 Id > m, Ree = 1/2}. 

(4.11) 

To lighten notation set r = ^{t) and suppress the t- and r-dependence from the 
notation. Let 9q = arcsin(l/2r) and notice 6o 'x 12 as r —>■ oo (t \ 0). For small 
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9 we have 


ire‘® -(7(re*®) =i/>(i)--r2g"(r)(e'® - if + - l)^) 


= ^(t)--r^\q''{r)\{e^ + 0{e^)) 


(4.12) 


The estimate is uniform in r = ^it) by Assumption |2.3[iii) . Let e{r) \ 0 with 
£{rfr‘^q''{r)/\ogr —^ oo (this is possible by Assumption |2.1[ ), then 

f exp(tre‘® — (/(r e*® ))ir e*® d0 ~ iy/27r|?//^(t)|e’^^*^ . (4-13) 

J-e(r) 2 


2 i 


As we veer away from r = ^{t) along Ti, the real part of ^t — q{^) decreases. Indeed 
for 0 G (0, 7r/2) 


—Re(tre‘® -^(re*®)) = Re - q {^))) 

Qf/ 

= -Im {^q (r) - ^q i^)) 


^—r 


< 0 . 


(4.14) 


At the very end we have used Assumption |2.3[ ii) It follows that 

^ f ° exp(tre‘® - g(r e*® ))ir e*® d0 
Js(r) 


rir 

< - exp 


1 


- + o 


= o(^ex]: 

Taking complex conjugates, we obtain a similar estimate for the integral from —9o 
to —£{r). Together with (4.13) we obtain 

1 /" eti-g(i) _ -y/2Tr\^"{t)\ . (4.15) 

2 i Jr, 2 


It remains to estimate the contribution from r 2 . Because of the monotonicity (4.14) 
we have 

Re(t^o - q{^o)) < W) - \T'‘^W'f)Hrff + o{l)) (4.16) 

Assumption |2.3[ii) ensures that 




d? 


< 


exp^Re [t^o - g(Co)] +0(logr))) = o(exp(i/;(t))) . (4.17) 

□ 


5. Critical point and Gaussian approximation 

Here we prove Lemmas |2.4p.7| and Proposition |2.12| and we address steps 4 and 5 
of the proof strategy. 

5.1. Variational analysis of fn{x)- Critical scales. 


Proof of Lemma [114 We treat the case a = 0. Under Assumption |2.1| (i), q' is 
strictly convex and decreasing. Therefore 

Nn-X 


fnix) = fix) - 


na^ 

j/ > 


>fiO + f'iOix-x*J- 


Nr,,-. 


na^ 


m - iv.. 




(5.1) 
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with equality if and only if a; = x*. If Nn < N*, we obtain /^(x) > 0 on (0,oo) 
and parts (a) and (b) of the lemma follow right away. 

If Nn> N*, then /^(x^) = {N* - Nn)/{na‘^) < 0 and lim^^^oo fni^) = oo, so by 
the intermediate value theorem has at least one zero in (x*,oo). On the other 
hand 


fnix) = q'\x) + (ncr^) ^ = g"(x) - g"(x*) > 0 on (x*,oo) 


(5.2) 


so /(j is strictly increasing and has exactly one zero x„ in (x*,oo), moreover 
fn(xn) = min(a;._oo) fn- Since g'(x„) > 0 by Assumptionand q'{xn) = (iV„ - 


x„)/(ncr^), we must have x„ < fV„. This proves the first part of (c). 

If in addition to Nn > N*, we have limsup„_^ 3 o Nn/{na^) < lima;^^^ q'{x), then 
limj-N^o fn{x) > 0. We have already observed that /(i(x*) < 0. By the intermediate 
value theorem, f'n has at least one zero x’n in (0, x*). Since /(((x) = q"(x)—q"{x*n) < 
0 on (0,x*), the zero is unique and corresponds to maximizer. This completes the 
proof of (c). The proof of (d) is similar to (c) and therefore omitted. □ 


Proof of Lemma [ITdt We treat the case a = 0. Write fn{x) = In[x,Nn) with 
In{x,y) = q{x) + [y - xf/[2na'^]. For y > Nf, let Xn{y) > x* > x'„(j/) be the 
solutions of dxln{x,y) = 0, with x'niy) well-dehned for < na^snpq' only. No¬ 
tice that X I—>■ I„{x,y) is increasing in (0,x(j(j/)), decreasing in (x(j(2/),x„(j/)), and 
increasing in (x„(j/),oo). We have 


^ [In (x„ (j/), y) - /(0,1/)] = -— 
ay '■ ■' na^ 


liy) 


na^ 


na^ 


< 0 . 


(5.3) 


As y \ N* at hxed n, a careful examination of the proof of Lemma |2.4| shows 
Xn{y) \ X*n and x'niy) /• Xniy), henCe IniXniy),y) IniXn.Nf). But X 
Inix,N*) is strictly increasing on (0, oo) because for iV„ = N*, d^Ini'^Nn) = 
f'nix) > 0 by Eq. (5.1), hence /„(x*,iV*) > J„(0, A^*) and by continuity 


hm [/„(x„( 2 /),y) -/(0,y)] > 0. 


Assumption 2.1 implies that q{y) = oiy) as y 


oo. It follows that 

,,2 


lim [/„(?/, y) -/(0,i/)] = lim [qiy)- 


y 


y—^oo 


2na‘^- 


= —oo. 


(5.4) 


(5.5) 


Eqs. (5.3)-(5.5) guarantee the existence and uniqueness of a solution y = N* 


to the equation /„(x„(i/),y) = /(0,y), and (a)-(c) follow with the observation 
fnixn) - /n(0) = [/„ (x„ (y) , ?/) - IniO, y)]\y=N^ ■ □ 


Proof of Lemma \2.6[ As noted after Assumption 


2.1 we have x"^q"[x) = 

—oo, moreover from the definition (2.10) of x* and the observation x* —>■ oo we get 

^2 


1 = lim na \q"iXn)\ » 


na 


(5.6) 


hence x* » The inequality x* < N* follows from the definition 2.10 of 


and the positivity of q'. The inequality Nf < Nf* holds true by dehnition of N* 


By Assumption 2.1 v) there exists C > 0 such that ^(x) < Cx°‘ for all sufficiently 
large x. Fix C > C and Nn > (2C"ncr2)i/(2-“). Then for large n. 


qjNn) 

Nn/i2na^) 


< 2na^CNf 


C , 


(5.7) 
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Write CjC = 1 — e. It follows that fn{Nn) = q{Nn) < (1 — ^ (1 ~ £)(1 + 

o(l))/„(a), and a fortiori min/„ < /„(iV„) < fn{a), which shows > N**. This 
proves N** = and completes the proof of the first part of the lemma. 

Next suppose by contradiction that N*lx. 
have na^q'{x^)jxl 


0. Since x* —>■ oo, Assumption 


1. T hen by Eq. (2.10) we must 
yields na^q''\x*^) 


2.1 


contradicting the definition of x* . Similarly, the assumption N*/x* 
na^q"{x^) —>■ oo, contradicting again the definition of x* . So N*/Xn stays bounded 


away from 1 and from oo. 


0 , 

oo leads to 
ed 
□ 


Proof of Lemma [1/^ Remember q'{x^) = [A^* — x*]/[ncr^] by definition of N*, and 
q'{xn) = [Nn — x„]/[n(T^] by definition of x„. Since q' is strictly decreasing we 
deduce 

N - X N* - X* N* 




= q{Xn) < q{x*J = 


< 






(5.8) 


so Nn — Xn < N* (see Figure]^. Since q' is strictly convex, we have q'{Xn) > 
q'{xn) + q"(x„)(x* - Xn) hence 


> 


q'{Xn) - q'[x*n) {Nn - Xn) “ (IV* “ <) 


0{NI) 


Xn-X*n na'^(Xn-X*) 

We also know that q"{xn) < 0, so we obtain 

1 1 


fn(.Xn) = q"{Xn) + 


na^{Nn + 0{N*n)) 


(-O’ 






(5.9) 

(5.10) 

□ 



(a) Nji < N* : no solutions 



(b) — N* : one solution 



(c) N-n > N* : two solutions 


Figure 3. Solutions to q'{x) = {Nn — x)/{na'^) (=critical points 


of fn) as in Lemma 2.7 


For the proof of Theorem |2.9| in the case Nn ~ N* we need the following. 


Lemma 5.1. Let W, = N* We have 


(Ar*)2 

for some e > 0 and all sufficiently large n, 


2na'^ 


(5.11) 


Proof. Assume a = 0 and g(0) = 0. By Lemma 2.5 we already know fn{Xn) — 


fn{0) > 0, we prove /„«) - /n(0) > sfn{0). Since /(.(x*) = 0 for Nn = N*, 


fu{x)= fn{y)dy=^ {na'^\q”{y)\-l)dy. 


na 


(5.12) 
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For 2 / > ai*/2 the integrand stays bounded away from zero, hence 

AW > 

n<7^ 

for all X > x*/2 and some c > 0. Then 

/„K)-/„(0) = 


fn{y)'^y > 


SncT^ 


(5.13) 


(5.14) 


and the statement follows from TV* = 0(ai*) (Lemma 2.6 1 and /„(0) = ^^"2 (1 
0 ( 1 )). The proof for a > 0 is based on a similar estimate of /n(a;*) — /n(2a) an 
therefore omitted. I 

Proof of Prop. \2.1^ Fix r S Nq. We have 


fLix) = /n(a^) 


(5.15) 


, Nn - X . 
na^ 

Clearly f^riNn) = q'iNn) > 0. Fix d G (0,e) and x G (x*, (1 + i5)a;*). Remember¬ 
ing ( 2 . 10 ) and the monotonicity of q', we obtain 

fUx) < - (1 + 0 ( 1 )) ^"-^^ + ^)^" 

" 'Yirr^ 


na^ 

1 

na 

1 


(TV„ + o(TV„) - K - 5x1) 


(5.16) 


<-^((l + o(l))£iV:-<5<) 


which is eventually negative because of cc* < TV* and 5 < e. It follows that f'^^ does 
indeed have a zero Xnr which lies between (l + i5)a:* and TV„. On (a;*, (l + 5)a;*) fnr 
is strictly decreasing by (5.16), on ((1 -f S)x'^,Nn) the second derivative satisfies 

^ .. 


fnrix) = -^(l-na‘^\q"{x)\ + 0 (—)] 
na V n / 


which stays bounded away from 0, hence fnr is strictly convex. It follows that Xr 
is the unique zero of f^r and the maximizer of fnr in {Xn,Nn). Moreover 


(l + n(l)) 






= q'{Xnr) < g'«) = 


TV*-, 




(5.18) 


hence Xnr = ^n — (1 + o{\))na^q'{Xnr) = Nn + 0(TV*). The asymptotic expression 
for fnrixnr) is easily checked. □ 

5.2. Critical points. We look for critical points { tn ,^ n ) with \ 0 and fn G 
(a, 00 ) of 

^ n ( t ,0 =- q ( 0 + n ( p ( t ) - mt + t ^ (5.19) 

(remember m = /xn+TV„). Since our contour integrals involve integrals over a = e*, 
it is convenient to work with functions of a single variable t: for t < limj,_>(i \q"{x)\^ 
let 

4'n(T) = = n(p{t) -mt + (5.20) 


where f{t) is the solution of q'{^{t)) = t as on p. 17 Then (t,^) is a critical point 
of <i>„(t, ^) if and only if ^ = ^{t) and 4'(j(t) = 0. So instead of looking for bivariate 
















22 


NICHOLAS M. ERCOLANI, SABINE JANSEN, DANIEL UELTSCHI 


critical points, we may look for zeros of in (0,sup |g"(a:)|). For later purpose 

we note the relations 


-m + = -{Nn - ^{t)) + na'^til + 0{t)) (5.21) 

= 710-2(1 + Olt)) + (5-22) 


with il}(t) defined in (4.9). The variable t and the analysis of are in some sense 
dual to the variable x (or and the variational problem fn{x) = min analyzed 
in Section 5.1 The analysis becomes more involved, however, because we need to 


take into account correction terms from J 2 j >3 KjP(j\. 


Lemma 5.2 (Inflection point of Let S > 0 such that inf(o_ 5 ) ip'^it) > 0. Then 
for all sufficiently large n, has an inflection point t* € (0,(5). Any inflection 
point satisfies ~ q'{x'^), and ifn is concave below the smallest inflection point 
and convex above the largest inflection point. 


If the inflection point is unique, then ikn is concave on (0,t*) and (t*, S). In general 
we do not know whether the inflection point is unique, however the asymptotic 
behavior ~ i® uniquely determined, and all statements below hold for 

every inflection point . 


Remark 4. It follows that t* is of the order of N*/n: from Lemma 5.2 and the 
definition of N* we have t* ^ 

n ■ n n./T^ 


and then Lemma 


2.6 


yields 


N* 

r n ^ j.* 


N* 
< C-^. 
n 


(5.23) 


Proof of Lemma \5.S\ As a preliminary observation, we note that any solution of 
'i'nit) = 0 converges to zero: this is because —' 0 "(t*) = n(p''{t) > ninf(o,j) >p'' —t oo. 

As t \ 0 at fixed n, we have ip^it) —)■ cr^ and fi'^t) —>■ —oo hence 4'"(t) —)• 
—oo. On the other hand we may choose e„ in such a way that e„ 0 and 
4'"(eji) —>• oo: indeed 'ij)"{q'{x*^)) = l/q"{xf) = —tkt^ by dehnition of a:*, so 
choosing S'^ch a way that |i/)"(c„)| ^ lif"{q'{x’^))\ = na^ we find 

= (1 + o(l))ncr 2 + o{na‘^) —>• oo. It follows from the intermediate value 
theorem that 'I'"(t) = 0 has a solution t* in (0,e„). It satisfies 


g"«) 


na^ 



— 1 + 0(t*) -t 1. 


(5.24) 


Assumption [2^ iv) and its consequence (2.5) imply that f,{Tf)/x*^ —>■ 1 and t* = 

q'mi)) - g'«). □ 


4'(j(t) is positive for small t and decreasing on ( 0 ,t*) and increasing on 
Define = C(C) notice ^ a;* from the proof of the previous lemma. 

For Nn <C n, let ijn Nn/{na‘^) be the solution of (3.4). Notice 4'(j(ry„) = 
00 . 


Lemma 5.3. Let 6 > 0 be as in Lemma 
ikn. Let Nn —t 00 with (1 + s)N* < Nn n. Then for sufficiently large n, ^'n 
has exactly two zeros in (0,(5), one zero tn G (0,t*) and another t'n € (tn,r]n)- Set 
f.n-=^{tn). We have 

0<Nn-^n = O{N:), 

limsupn^oo tn/tn < 1; 1™inf„^oo > 1, and liminf ^/? 7 „ > 0. 


5.2 s > 0, and tn an inflection point of 
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Remark 5. When Nn ^ N*, we use Eq. (2.4), Lemmas 2.6 and 5.2 and find 

tn = q\^n) - « g'«) = o(^) 




(5.25) 


hence = o(i*). When iV„ = 0{N*), we have instead < ct* for some 

c > 0 and all sufficiently large n: The upper bound is part of Lemma |5.3| For 
the lower bound, we note that = q'{^n) > q'{Nn) because q' is decreasing and 
in < Nn- Since is of the order of A^*, Eq. (2.4) shows that q'{Nn) is of the 
order of q'{N*) which in turn is of the order of q'{Xn) ~ tj)- 

We check first that t* < r]n- With C > 1 as in Lemma 2.6 

n:-x. 


Proof of Lemma [X 
we have 


C ^ ?'«) = 




<{i-c-^) 


N* 


na^ 


so 


1-C 


-1 


1 + e 

• oo, and as n 


(5.26) 

(5.27) 


(5.28) 


tlhn < (1 + o(l))(l - C-^)N*JNn < (1 + 0(1)) 

is bounded away from 1. As t —)■ 0 at fixed n, dt^t) ^ i{t) - 
^niVn) —t oo. Furthermore 

K{t*n) = -{Nn - C) + (1 + 0{l))na\'{i*n) 

= -{Nn-{l + o(l))0 + (1 + o(l))ncr^g'«) 

= -{Nn - (1 + o(l))<) + (1 + 0 (l))(iv: - O 
=-(7V„ - TV*) + o(7V*) ^-00. 

The intermediate value theorem proves the existence of a zero G (0;^n) 
another zero t'n S {tn,r]n)- 

Suppose by contradiction that tn/C —>■ 1- Then the identities = <f{in), 
tn = q'{in) and Eq. ( |2.4[ ) imply in ^ in ^ x* and an estimate analogous to (5.28) 
shows 4>(j(t„) —)• —oo, in contradiction with 4>(j(t„) = 0. It follows that t„/f* and 
in/Xn are bounded away from 1. In addition, tn <t*^ = 0{N*/n) and 

Nn-in^naHn = 0{N*). (5.29) 

For the lower bound of t'n, we use t'n > tn q'{^V) and Assumption 2.1 to get 
in < (1 + o (l)X where q'{i'n) = t'n- Since ^ = 0(??„) ^ 0, 'i’n{tn) = 0 together 
with (5.21[) yields 


t 


Nn-i'n 


71 ( 7 ^ 


{'-k)M'-^) 


(5.30) 


and A^„ > (1 + e)N* implies liminf> 0. Notice that, in view of ( 5.27| ), we 
have a fortiori liminf > 0 . 

We have actually shown that for every inflection point , t/t* < 1 stays bounded 
away from 1. In particular, in case of non-uniqueness of t* we may choose 
as the smallest inflection point of 'i'n- Then i>'^ is strictly increasing on (0,t*), 
consequently the zero is unique. A similar argument shows that t^ is unique. □ 

When liminf„_>oo Nn/n > 0, the sequence ijn is either no longer defined or it does 
not converge to zero. The previous lemma is modified as follows. 


Lemma 5.4. Assume Nn —>■ oo with liminf A^„/n > 0. Then there exists Jq > 0 
such that for large n, itn has exactly one critical point tn in (0,(5o)- The critical 
point lies in (0,t*) and it satisfies 0 < AT„ — i{tn) = 0{N*). 
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Proof. The existence and uniqueness of a critical point in (0,t*) as well as the 
properties of 5(t„) are proven as in the previous lemma. Fix Jq > 0 such that 
p'{So) < /J + 5 liminf(7V„/n) =: fj, + e/2. Then 

- 'F'(^o)) + C('^o) < -«(| + o(l)) + C(^o) -oo- (5.31) 

It follows that < 0 on (t*,(5o). □ 


5.3. Hessians. Let G (0, t*) be the critical point of 'l'n(t), and = ^{tn)- Thus 
(trij^n) is a critical point of Lemma B.l shows 

det Hess (tn, Cn) 1 - nRe ) | q” ) | 


■ W'iUl 

Lemma 5.5. Assume Nn —>■ oo with liminf„_>oo 


(Nn/N*) > 1 and let be 


the unique critical point o/$„ in (0 ,t*) x (a, oo). 

(a) If Nn » N*, then det Hess ^n) —)■ —1. 

(b) If Nn = 0{N*), then detHess$„(t„,^ri) =-(1 - n(T2|g"(^„)|) + o(l) and 
it stays bounded away from zero. 


Proof, (a) If Nn ^ N*, then by Lemmas 5.3 and 5.4 we have fn = Nn + 0(IV/) ^ 
Nn, hence in particular 3> N* > a;*. Exploiting the monotonicity and the 
convexity of q', we have 

0 > 9"(C ) > ~ 

from which we get 


/N*\ 


1 

ncr^ 


fn-K A^„(l + o(l)) 

det Hess $„(t„,^„) = -1 + ncr^(l + 0(t„))|g"(^„)| = + 


(5.33) 


-1. (5.34) 


(b) If Nn = 0{Nn): By Lemma 5.3 > (1 + s)Xn for some e > 0. Consequently 

|9”(^n)| < (1 ~ ^)\q"{^n)\ = (1 ~ ^)/(^o-^) for some 5 > 0 and large n, from which 
we deduce that 

1 - nReq>''{tn)\q"{^n)\ >!-(! + 0{tn)){l -5) = 1-6 + 0{tn), (5.35) 

in particular the expression stays bounded away from zero. We also have 

det Hess $„(t„,C«) =- na'^\q"{^n)\) + o{tn)n(T'^\q"{^n)\ (5.36) 


Since > x*, the estimate (5.33) still holds true and 

'Nf 


tnncj'^\q"{^n)\ = = 0{tn) 0 . 


(5.37) 

□ 


5.4. Gaussian approximation. Next we address the Gaussian approximation for 
the evaluation of i7„. Because of Theorem |4.4[ we need not deal with a bivariate 
integral and instead may use 


1 

TT 


gnRov(i)-mi IjnG(e*)dG 




V'|i/’"(t)|e'^"(*^ dt (5.38) 


as n —>■ oo and \ 0. Remember that 'I'„(t„) = ^n{tn,^n) and from Eq. (5.32) 
and Lemma 15. 5 1 


i’ujtn) ^ 1 

4'/(t„) l-na‘^\q”{^n)\ 


(5.39) 
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with a denominator bounded away from zero. The following technical lemma helps 
estimate the prefactor (t)\. Set 


Rn{t) := 'i’nit) + log \/|i/i"(t)|. 


(5.40) 


Lemma 5.6. Let N„ 


oo with liminf(7Vji/A^*) > 1 and G the zero of 

from Lemma 5.3 Then R'„ has at least one zero Sn G (0,<„). Moreover there 
exists a sequence \ 0 such that every such zero lies in ((1 — 


Proof. As t \ 0 at fixed n, using Lemma [4.3[ we have 
1 


<W = <(i) + A- 


= n(Re y}'{t) - /i) - iV„ + f;{t) + ^ 


= n{a'^t + Ofi^)) - Nn+ ip'{t) + 0{l/t) 
= -iV„ + o(l) + (1 + o(l))V''(t) -t oo, 
hence 1 1 —>■ is initially increasing. At t = we have 






< 0 . 


(5.41) 


(5.42) 


2 fj'ftn) 2 i/>"(t„) 

The intermediate value theorem guarantees the existence of a zero s„ of i?(j. Set 

1 1p^^\sn) 


yn ■■= -fl'iSn) + 


2 V'"(s„) 


(5.43) 


In view of Lemma 4.3 we have yn ~ f^'isn) and by Eqs. (2.4) and (4.10), q'{yn) 


Sn. From Sn < tn < t’^ and Lemma [27f| we get that is larger than x* and actually 
bounded away from it. By the definition of s„, 

0 = n(Re(/?'(s„) - y) - Nn + yn = na‘^Sn{l + 0(s„)) - [Nn - yn 

hence 

Un 




q'ivn)- 


(5.44) 

(5.45) 


Let Xn € {Xn, Nn) be the solution of q'(xn) = {Nn — a:„)/(ntT^). Now 
d / ,, , Nn-y\ _ ^ , 1 I - na'^\q" {y)\ 


dy 


{q'{y) 


= q"{y) + — = 


(5.46) 


and na’^lq"{y)\ stays bounded away from 1 when y > (1 + S)Xn, we find that for 
some c > 0, we have 


Q'{yn) - 


Nn-y 


> 


c\yn - Xn 


(5.47) 


hence yn — Xn = o{Nn — yn) = o{Nn). If fV„ » TV*, Lemma 2.7 says Xn Nn and 
we deduce yn — Xn = o(xn). If TV„ = 0{N*), we use x„ > x* in conjunction with 
Lemma 2.6 and find yn — Xn = o{N*) = o(x*) = o(x„). 

Thus we have checked that i/„ ~ x„, which in turn yields Sn q'{xn). An 
entirely similar argument shows ~ q'{xn), so we must have Sn ~ In- This holds 
for every zero in (0,t„), in particular the smallest one, and the lemma follows. □ 


Lemma 5.7. Let Nn 

rt' 


OO with liminf(A„/A^*) > 1 and Nn = 0{N*). Then 


\/|V'"(I)| dt' 


2'k 


1 - na^q'f^n) 
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Proof. By Lemma 5.2 and |5.3[ is increasing on (0, and decreasing on (i„, 

We use a Gaussian approximation to around and adapt m Lemma 2.1]. 

First we check that the window (t„ — +en) contributing most to the Gaussian 

integral fits amply into i.e., 

= o(min(t„,4 - in))- (5.48) 


<^n • — 




By Lemma 5.3 we have < (1 — ^)i* < (1 — for some S G (0,1), hence 


t'n~tnP Stn and 


min(<„,t(, - tn) > 


(5.49) 


Eq. ( 5.39| ) and Lemma 5.5 sho w tha t 'I'"(t„) is of the order of i/>"(t„), of the 
order of l/|t/’n(in)|- By Lemma 4.3 i^|V'n(in)| —t oo hence t 0 and (5.48) 

follows. The same argument shows that for 

C„ —>■ OO with CnSn = oftn) (5.50) 

we still get c„e„ = o(min(t„, — in)) ■ Second, we observe that the prefactor tp”{t) 
is essentially constant on the relevant window: Because of (u)/tp"(u) = 0{l/u) 
(Lemma |4.3[ ), we have 

rit) 


V'"(in) 




and 


supl 


rit) 


- 1 


|i - ini < C„£„| = 0(c„£„) 0. 


I V'"(in) 

Third, we note that cubic corrections can be neglected: 

4'n(i) = ^nitn) + ^ ^(1)) ^n(in) (i — tn)"^ 

uniformly in |t — t„| < c„£„. To this aim write, with the help of (5.52), 
nt) = niritn) + 0(c„£„)) + ritn)i^ + 0(c„£„)) 

= <(!„) +C(in)(o(^gg) +0(c„£„)) 

Now nlf}"(tn) = nq"{^n) = 0(1) by Lemma 
Eq. (5.39|) and the same lemma, hence 


(5.51) 

(5.52) 

(5.53) 


5.5 


and V’"(i„) = 0(4'"(t„)) by 
Kit) = {1 + OicnsrKitn) (5.54) 


in |i —i„| < c„£„ and (5.53) follows. Eqs. (5.52), (5.53) and (5.39) yield 

ptrL-\-Cnen 


^\rit)\dt - 


' tn—CnSn 


l2Trritn] 

Kitn) 


27r 


1 - na‘^q"ir) 




(5.55) 


Our next task is to estimate the integral on (0, tn — CnSn) and (in +Cnen, in), taking 
into account that the prefactor \/\'>p"it)\ goes to infinity. On both intervals we have 


4'n(i) < ^'n(in) “ 2^^ 0(c„£„))c^. 


(5.56) 


For i > i„ + c„£„ we have —ip''it) < —V'"(in) hence 

\/|V'"(i)|exp(Tn(i)) < exp(^'„(t„) - ^{1 + oil))cl + ^ log |V'"(in)|) • (5.57) 
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By Lemma |5.6| we may choose c„e„ 0 in such a way that is increasing on 
(0,t„(l - c„e„)), which yields 


+ log \/|V'"(t)| < - CnSn) + log “ C„e„)|. 


(5.58) 


This estimate, combined with Eqs. (5.52) and (5.56), shows that (5.57) holds true 
not only in (t„ + but also in (0,t„ — c„e„). Next we check that we can 

choose c„ —>■ oo so that not only c„e„ = o(t„) but in fact 


- (1 + 0(1))C^ + log -OO. 


(5.59) 


By Lemma 4.3 (iii), as t \ 0, we may estimate \'4’"(t)\ as follows: fix tp > 0 and 

rto 


take t G (0, to)j then 

log(-i/>"(t)) = log(-'i/>"(to)) - 


rio 


to 


du < logi-ip (to)) 


C 


-du 


= log(-V'"(to)) + Clog y = 0(1 logt|). 

In particular log |i/)"(t„)| = 0(logt„). On the other hand 

Cn£n\^Kitn) {t„)\ 


I log t„ 




|l0gt„ 


(5.60) 


(5.61) 


(5.62) 


The ratio 4'"(t„)/i/!"(t„) stays bounded away from 0 and by Lemma 4.3 

—,-,- > OO. 

|logt| 

Thus we may find a function a;(t) —0 such that a;(t)^t^|'!/'"(t)|/| logt| still goes 
to infinity as t \ 0, set c„£„ = t„a;(t„), and then Eq. (5.59) holds true. The 
bound (5.57) then shows 

= o(e'*'’'(*-)). (5.63) 

Jo Jt„+c„er, ^ ' 


Eqs. (5.55) and (5.63) complete the proof of the lemma. 


□ 


Lemma 5.8. Let Nn 


OO with N* <C Nn <C n. Then 


rt' 




Proof. The proof of Lemma 5.7 applies without any changes, the end result simpli¬ 
fies because 1 — na’^q"{f^n) —>■ 1 by Lemma 5.5 


□ 


Lemma 5.9. Let Nn —>■ oo with liminf(7V„/n) > 0. Let Sq > 0 as in Lemma 5.4 
Then 

rSo 


( \/|V'"(t)| dt ^ _ 

Jo 


Proof. By Lemma 5.4 4>„(t) has a unique critical point tn in (0,(5o) and we may 
use a Gaussian approximation on this interval—there is no need to restrict to an 
interval whose length goes to zero. Apart from this difference, the proof is identical 
to Lemmas [5^ and 15^ □ 
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6. Evaluation of contour integrals. Proof of the main theorems 
The proof of the main theorems starts from the decomposition 

F{Sn=n^l + Nn) = Hr,+Vn ( 6 . 1 ) 


where and 14, are defined as in ( |3.7| ). The correctness of (6.1) is checked as in 
Section [3 buil ding on the properties of G{z) proven in Section |4j For the proof of 
Theorem 2.10 it is convenient to decompose Hn further as Hn = where 


Hl = - 


gnRev(i)-(Mn+N„)i sin(nlm V9(f))dt 


( 6 . 2 ) 


and is a similar integral, but with integration from t'„ to rjn- In the proof 
of Theorem |2.11 we adopt slightly modified definitions and replace the sequence 
r]n ~ Nn/iria'^) in the domain of integration by another sequence \ 0 or by 
some fixed small e > 0. 

6.1. Evaluation of 14. Theorem 2.11| only needs upper bounds for 14, provided in 
Lemmas |6.1| and |6.2| Theorems |2.9| and |2.10| requires the full asymptotic behavior 
of 14 proven in Lemma |6.3| 


Lemma 6.1. Suppose A4 —t oo and Nn > Sn for some 5 > 0 and all n G N. Then 
for suitable C = Cs > 0, every sufficiently small p > 0, and all n gN 


1 

TT 


^nRe(p(ri+ie)-(fm+N„)ri ^0 ^ Q-Cnr] 


Proof. Let 5+ be the half-strip {t S C | Ret > 0, Imt G [0, tt)}. By Theorem 4.2 
we know that as t —>■ 0 in 5+ 

Re {ip{t) — fit — ^^t) = —^^Ret + ^cr^((Ret)^ — (Imt)^) -I- 0{t^). (6.3) 

Then for sufficiently small Ci > 0 and all t G S+ with max(|Ret|, |Imt|) < ei, the 
right side of Eq. (6.3) is smaller than —SRet/2, which shows 


1 

TT 


^nRe-fi{7]+W)-{(tn+N„)r] 0^0 ^ ^ ^-nSr]l2 


(6.4) 


for all T] G (0,ei). On the other hand on the unit circle G(z) is given by a power 
series with strictly positive coefhcients and therefore |G( 2 ;)| has a unique maximum 
at z = 1. It follows that for all t = iff with 0 G [ei, 27r—ei], we know Re (p(t) < 0. By 
continuity this extends to some thin vertical strip Imt G [ci, 27r — £i]. Ret G [0,£2] 
so that Re((^(t) — fit) < —fiRet and 


1 

TT 


^nRo(p(r/+ie)-(fj.n+N„)r] ^0 ^ 


TT — £i _ 


nfirj 


for all rf G (0,£2). To conclude, we let Gs ■= min((5/2,^). 

Lemma 6.2. Assume Nn = m — nfi -G oo and Nn = o(n). Let £„ 
= (1 + o(l))^^- Then for suitable constant G > 0, as n ^ oo, 


(6.5) 

□ 

0 with 


1 

/ ^nRetp{er,+iS)-me„) ^0 

271 Jr, 


< g-(l-So(l))Ar„e„/2 _ 
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Proof. The proof is analogous to Lemma 6.1 We start from the estimate (6.3). In 


a sufficiently small ei-neighborhood of the corner 0 of the half-strip 5+, we have 

Re — —t) < — icr^Ret^Ret. (6.6) 

Notice that El can be chosen n-independent: we only need—CT^(Imt)^+0((Imf)^) < 
0 for |Imt| < El- When Ret = £„ with £„ < the upper bound in (6.6) is in 
turn bounded by —(1 + o(l))fV„e„/2. 

When Ret = £„ is small but Imt is bounded away from 0 and 2tt, we estimate 

Re - ”i) < -^Rei = -(p + (6.7) 


n 


n 


n 


n 


and we conclude as in Lemma l6.1 


□ 


The proof of the next lemma is closely related to the treatment of moderate 
deviations for random variables with generating functions analytic beyond z = 1 
given by Ibragimov and Linnik [12) . 

Lemma 6.3. Let > oo along ^/n <C = 0(ri}~^) for some 7 > 0. Define 
Nn/{na'^) by (3.4) and I 4 , as in (3.7). Then Eq. (3.11) holds true. 


n 


Proof. Since Reip'{t) = fJ, + cr^t + O(t^) is strictly increasing for small t > 0 and 
Nn/n —0, we may fix <5 > 0 small enough so that for large n > ns, the equa¬ 
tion (3.4) has indeed a unique solution rjn € (0,5), which satisfies ~ Nn/{na'^). 


Arguments analogous to the proof of Eq. ( |6.5[ ) show 

1 cos(nIm7)(77„-ki6»))d6l 

Js 

We also have, uniformly in s e (0,5„), 

PiVn + is) = <p(?7n) + </5'(?7n)is - ^<p"(?7n)s^ -|- O(s^) 




( 6 . 8 ) 


= PiVn.) + + ^)s - Imv 3 '( 77 „)s - + O(s^) 


(6.9) 


and therefore 


Re <p(? 7 „ -I- is) - ^/i-l - = Reip{rin)- -if” {qn)s^ {I + 0{s))-Imip'{qn)s. 

( 6 . 10 ) 


Since ip"{qn) cr^, standard arguments show 
^5 

gnRe(//(7)„)-f Re(73"(77„)s^(l-|-0(s)) ^ 


1 

TT 


'0 


exp(nReyi(7?„)) 

V2TTna^ 


( 6 . 11 ) 


moreover the contribution to the interval from s > 5„ := (logn)/Y^ is negligible 
and Eq. ( 6.11| ) holds true with S replaced by 5„. By Theorem 4.2 and the relation 
ImG'(e‘) = |G( e‘)|Im 7 )(t) the imaginary part of p'{qn) vanish faster than any 
power of qn = 0(n~^), hence nhnip'[qn) can be neglected; the same argument 
works for nhnp"{qn). For the cosine, we look separately at (0,5„) and (5„,5). On 
(0, 6n), again by Theorem 4.2 supsg(o, 5 „) |Ini< 7 ’(ry„ -I- is)|, vanishes faster than any 


power of max(?7„,5„) hence it can be neglected. On (5„,5) we simply bound the 
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cosine by 1. As the contribution from to the integral (6.11) is negligible, 

combining with (6.8) we find in the end 

TA , ,,,,exp(nRev3(?7„)-^n-A^„) 

= (1 + o(l))— 


V 27171(7^ 

By Theorem |4.2| and Definition! 2.8 [ we have 


/ N \ 
nRe(^(?7„) - ^xn- Nn= (^1 + 0(^)j 


+ 0(e-”'^’'"). 

(6.12) 




(6.13) 


with correction terms expressed in terms of the Cramer series. In particular, the 
exponent goes to —oo as —N'^/n, i.e., slower than the term —n/i77„ = —Nn^x in the 
second term. Therefore the second term in Eq. (6.12) is negligible compared to the 
first and Eqs. (3.10) and ( |3.11| ) hold true. □ 

6.2. Evaluation of Here we focus on the regime liminf Nn/N* > 1; the case 
limsuplV„/A^* < 1 is treated in the proof of Theorem 2.9 In order to apply the 
Gaussian approximation from Section |5.4[ we need to drop the sine and replace 
Im(/5(t) with G(e‘). 


Lemma 6.4. Let £„ \ 0 faster than some power of n, i.e., ne^ —>■ 0 for some 
p > 0. Then 

sin(nlm(/5(t)) ~ nImG(e*) 

uniformly for t G [0, e„]. 


Proof. We have 

ImG(e*) =Im = e^® Im (^(t) = |G( e‘)| Im 
We know that 

sup |ReG(e*) — 1| = 0{en) —>■ 0, 

t6[0.e„] 

and the imaginary part vanishes faster than any power, in particular 

sup |ImG(e* )| = 0(e5^) —>■ 0. 

te[0,e„] 

It follows that a similar bound applies to lvcitf{t). As a consequence 


(6.14) 

(6.15) 

(6.16) 


sup 

te[o,£„] 


nlTo.ip(t) — sin(nlmi^(t)) 
sin(nlm ‘p{t)) 



= O 



-G 0 , 


(6.17) 


Thus sin(nlm(/5(t)) ^ nlnupif), uniformly in [0,e„]. 
Ivnipit) ~ ImG(e*) uniformly in [0,£„]. 


Eq. (6.14) in turn shows 

□ 


The dominant contribution to the Gaussian integral in Lemmas |5.7||5.9| comes 
from windows of width o(t„) around by Remarkj^and Lemma 2.6 = 0(t*) = 

0{N*/n) = This latter bound vanishes like some negative power 

of n, hence Lemma |6.4| is applicable on the interval contributing most to the 
Gaussian integrals. Outside we use the inequality | sin(nlm(/3(t))| < nlrmplt) = 
n(l + o(l))ImG(e*), and we find: Eor Nn -G oo with liminf„_>oo Nn/N* > 1, we 
have 


Hi 






na'^q' 




(6.18) 
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For Nn —>■ oo with iV„ ^ N*, 


Hi 


ne 




(6.19) 


Hi is estimated in the proof of Theorem 2.10 and is not needed in the proof of The¬ 
orem 2.11 Finally for liminf Nn/n > 0 and i5o > 0 small enough as in Lemma 5.4 

<5o 


■ f ° sin(nlm(p(t))dt - . 

Jo 


( 6 . 20 ) 


Remember that ^'„(t„) = $„(tri,^n) = -fnr{xnr) + o(l ) and 1 - ncr^|g"(fn)| ^ 
1 — ncj'^\ q"{xnr)\ by the definition of and Eq. (2.22), so the right-hand sides 
in Eqs. ( 6.18| )-(6.20) correspond to the relevant contribution in Theorems 2.10 
and 12 . in 


6.3. Critical scale: proof of Theorem 2.10 Let oo with liminf iV„/iV* > 

1 and Nn = 0(n^/(^““^). Let ? 7 „ ~ Nn/{na'^) be the solution of (3.4), define iL„ 
and Vn as in d^ , and Hi and H^ as in (|6.2[) Thus we have 


lP('S'ra — n/r -|- Nn) — Vn + HI + HI- 


( 6 . 21 ) 


By Lemma 6.3 the asymptotics of is given in terms of the Cramer series as in 
Eq. (3.11). Hi is evaluated with Lemma 5.7 and Lemma 5.8 as the right-hand side 


of Eq. (3.22). The proof is complete once we show = o(14). 


On {t'n,r]n) the function 4'„(t) is increasing by Lemma 5.3 


sup '^n{t) <^nir]n) = (nReip{r]n) - infi +Nn)rin')+iljir]n). (6.22) 

O' r, 1 ^ r 




The term in big parentheses can be reexpressed with the Cramer series and is 
exactly equal to the exponent in the evaluation of Vn (see Eqs. (3.10) and ( 3.11| )), 
while tl{r]n) —oo. The prefactor satisfies 

sup ^ v/|?/’"(t)| < \/\i^''{t'n)\ = exp(^0(logt(j)) = e^p(o{logn)y (6.23) 


(remember (5.4) and tl > t*m with of the order of Nl/n 1/y/n). On the 
other hand rjn faster than some power of n, so by Lemma |4.3| we have for some 
constant C 

\tjjir]n)\ > |log? 7 „| > Clogn » log y|'0"(t'„)|, 


whence 


tp{'nn)+Vgn+ sup log-v/|V'^^(<y|-oo. 


In view of (6.22), (3.10) and (3.11), we obtain 

rv 


Hi = - 


e’"^®^(*)-("'"+^’')*ImG(e‘)dt = o(E„). 


which concludes the proof. 


(6.24) 

(6.25) 

(6.26) 
□ 
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6.4. Big jump: proof of Theorem 


2.11 


Let Nn —>■ oo with iV„ :g> 


Notice that, by Lemma 2.6 we then have Nn ^ N**. We distinguish two cases. 

Case 1: liminf„_).oo Nn/n > 0. Fix Jq > 0 as in Lemma 5.4 and define Hn and 
Vn as in (3.7) but with Jq instead of rjn- The asymptotic behavior of Hn is given 
by Eq. (6.20). The proof is complete once we check Vn = o(i7„). 


Vn is exponentially small in n by Lemma 6.1 Remembering Eq. (2.22) and (??) 
we get 

4'n(tn) = -fnr{Xnr) + o(l) = -q{Nn) + ^(1 + o{l))na'^ q'{Nnf + o(l) 

>-g(iV„) + o(l)>-C7V“. (6.27) 

It follows that Hn ~ nexp(4>„(t„)) goes to zero slower than exp(—cn“) for some 
c > 0, hence 14, = o(i7„) and 


P( 5 '„ =nn + Nn)'^ Hn^n , (5,28) 

If Nn ^ n, set e„ = t'n- By Remark]^ the critical point t'n G {ttiyVn) 

is bounded from below by some constant times rjn ^ iV„/(ntT^). We define Hn and 
14 as in (3.7) but with £„ instead of rjn- Hn is evaluated as in (6.19), which yields 
Hn ~ nexp(4'„(4)); Eq. (6.27) stays valid. I 4 is estimated by Lemma 6.2 which 
yields I4 = 0(exp(—A^„t(j/2))). Now Nnt'n 00 much faster than 1V“. Indeed t'n 
is bounded from below by some constant times Nn/n, hence 

N 


Nnt'n 


= 0{nN^-^) 


0 . 


(6.29) 

□ 


It follows that 14 = o{Hn) and Eq. ( |6.28[ ) stays true. 

6.5. Small steps: proof of Theorem |2.9[ The proof of Theorem |2.9| requires two 
more technical lemmas, proven at the end of this section. Remember the function 
Rn(t) from (5.40[). 


Lemma 6.5. If limsupn^oo 
large n. 


Nn/N* < 1, then R'n > 0 on (0,77„) for all sufficiently 


Lemma 6.6. Let Nn —>■ 00 with Nn ~ N*. Set f/{x) := q{x) + ■ Suppose 

that there are infinitely many n € N for which the equation R'nft) = 0 has a solution 
Sn G (0,??n)- Then Sn ~ tn and 

RniSn) = -/:«) + + 0 {\ 0 gn). 

The zero need not be unique—in case of non-uniqueness the lemma applies to 
every choice of s„. 

Proof of Theorem \2.9\ Let — >■ 00 along Nn < (1 + o{l))Nn. Let ri„ 

Nn/{na‘^) be the solution of (3.4) and define I 4 and Hn as in (3.7). By Lemma 


the asymptotic behavior of 14 is given by Eqs. (3.10) and 3.11 
verify that = o(I4)- We estimate 


6.3 


so it remains to 


ntjn 

tG(0,r)„) 


Hn<— sup sup e'^-W. (6.30) 

te{0,ri„) 


Just as in Lemma 5.6 one checks that Rn{t) is increasing for small t. We distinguish 
two cases. 
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Case 1: lim inf Nn/N* < 1. Then Lemma 6.5 shows, for large n, 

sup Rnit) < RniVn) = {nReip{r]n) - {nfl + N„)r]n) + tllitju) + log \/|'(/'"(?7n)l 
■.(r\ \ V / 


ie(0,77n) 


(6.31) 


and we deduce from (3.10) 
Hn ^ nr]. 

Vn. ~ TT 


(6.32) 


V2Trna‘^ ejipy-ip(rin) + log a/|i/’"(? 7 ™)| 

= exp(|V'(?7n) + log + O(logn)) 

which goes to zero by an estimate analogous to ( 6.24| ). 

Case 2: ~ N*. If R'^ reaches its maximum at t = the estimate (6.32) still 

applies. If along some subsequence {nj), R'„ reaches its maximum at some interior 


point Sn G (0, rjn), then we must have R'^{sn) = 0 and by Lemmas 6.6 and 5.1 

*2 


N* 

sup i?„(t) <-(l + e + o(l))-^+0(logn). 


te(0.r)„) 


2na^ 


(6.33) 


for some e > 0 and all large n. Since 14 = exp(—(1 + o(l))^^^ + O(logn)) we get 

(6.34) 




< exp(^-(e + o(l))^^ +0(logn)y 


By Assumption |2 .1 1 and Lemma [2^ 

1 I *11 ^ logn + 0(l) 

- 2 = ^ i^n)\ > -^ > 


na^ 


x: 


2 — 


Xy, 


(6.35) 


so logn = o(iV* /n) and the right-hand of (6.34) goes to zero. 


We have checked in both cases that = o(14), which concludes the proof. □ 


Proof of Lemma \ 6. 5| Suppose that the equation R'^[sn) = 0 has a solution s„ S 
(0, rin) for infinitely many n. Thus (s„) may be defined only along some subsequence 
(rij), which we supp ress from the notation. Define !/„ ^ V'^(sri) as in (5.43). Pro¬ 
ceeding as in Lemma 5.6 we find that q'{yn) ~ [Nn — Vn )/By the convexity 
of q' and the definition of iV*, we have 

(i-ko(i))- ^ = q{yn)> -^ =^^ (6.36) 

hence N* — Nn < o{Nn — yn) = o{Nn) and limsup iV*/A^„ < 1 i.e. liminf /iV* > 

1 . 

So if limsup„_,,go Nn/N* < 1, we must have R/ 0 on {0,r]n) except possibly 
for finitely many n. From the proof of Lemma 5.6 we know that lim^x^o R'n(f) = oo 
for all n S N, and Lemma [67^ follows. □ 


Proof of Lemma 6.6 For t G (0, rjn) with rjn ~ Nn/(na-^) we have 
R'n(t) = n(Re(p'(t) - jr) - Nn + 'if'{t) + 

= nah - (K - m) + o{nt) + {N/ - 1V„) + oQ 

= nah-{N:-m)+o{N:)+om)- 


(6.37) 
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In terms of the dual variable ^ = ^{t) = the equation = 0 reads 

na^'iO - (K -0 = oiK) + o(a (6.38) 


At ^ = X*, the left-hand side of (6.38) vanishes by definition of N*. It follows that 

ncr'^q'iO - {N*-^) = f {l + na'^q”{x))dx. (6.39) 

J tc* 

Fix 7 > 0. On ((I -|- 7 ) 0 ;* , 00 ), we have I — na‘^\q"{x)\ > I — na^\q"{{l + 7 )a:*)| =: 
c~f > 0 and for all x > x*^, the integrand is non-negative. As a consequence, 

na'^q'i.S.) - {N* - C) > - (I + qX) (6.40) 

for all ^ > (I-|- 7 )a;*. Suppose that Eq. ( 6.38[ ) has a solution with > (l-|- 27 )x* 
along some subsequence. Then N* = 0(a:*) = and 


C77C < c^iC - (1 + 7)0 = o{C), 


(6.41) 


which is a contradiction (remember > x* —>■ 00 ). It follows that for every 7 > 0, 
there are at most finitely many n for which > (l-|- 27 )x* , hence limsup^(j/x* < 1. 

The case ip'{r]n) < Cn < (1 ~ 27 )^* is treated in an analogous fashion, based on 
two observations: first, na‘^q''{x) -|- 1 < —c-y < 0 for all x S (a, (I — 7 )x*) and some 
Cj > 0. Second, since is of the order of N*/n i.e., of the order of t* ~ q'(x*^), 
the estimate (2.4) shows that is bounded from below by some constant times 
iV*, i.e., we still have N* = and —)■ 00 . We find liminf^(j/x* > I, hence 

altogether lim A/x* = I. This applies in particular to £' := tb'isn) i.e. s„ = q'(£,'„)■ 
Eq.pland tl ~ < 7 '(x;) (Lem ma yield lims„/f; = 1. 


By Remark and Lemma 
logs„ = O(logn) and by (|5.4[ 


2.6 


we have Sn. ~ > const iV*/n ^ 1/v^ hence 

(6.42) 


log|'!/'"(sn)| = 0 (logs„) = 0 (logn). 
Furthermore for t G (0,r/n) 

4'«(i) = ^n{aH^ + o{r]l)) - N^t + t({t) - q{^{t) 


{-9(C(i)) + - (X - + 0{{N* - A„)? 7 „) -h o{nril) 


(6.43) 


The two remainders are o{{N*)'^/n) by our choice of fV„. Write gn{t) for the term in 
curly braces. At t = q'(x*) = [N*-x*„]'^/[na'^] we have ^(t) = x*, g„(t) = -/*(x*) 
and (?(j(t) = 0. Moreover for t £ (0, ? 7 „) 

g”{t) = na'^ + C'(0 = na'^ - = 0{n(T^). (6.44) 

Here we have used that t < gn implies that ^(t) is bounded from below by a constant 
times N* or equivalently, x* and therefore \q''{^{t))\ is bounded from below by some 
constant times l/{na^). We deduce 


gn{sn) + fn{Xn)\ < ^(t„ “ (x*))^ sup |g"| = o{tlna^) = o{ngl) 

2 0 ,t 7 „ \ n J 


The estimate on Rn{sn) follows from Eqs. (6.42), (6.43) and (6.45). 


(6.45) 
□ 
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7. Extensions 

The methods described in this paper depend on being able to extend weight 
functions defined on integers to analytic functions on some region in the complex 
plane ( Assumption |2.2] ) . Here we explain how one might proceed when the existence 
of an analytic extension is not a priori given. 


7.1. Probability. Log-convexity. A first option, for weights satisfying Assump¬ 
tion 2.1 is to proceed as in m At the end of the section we comment on the role 
of log-convexity and a possibly broader perspective. As shown in m Lemma 2.1], 
Assumption |2 .1 1 is enough to guarantee that for every d > 0, as n —>■ oo, 

rS 


p{n) ^ [ e -\/27r|'(/)"(t)| dt, 

Jo 


(7.1) 


a statement in some sense inverse to Theorem |4.4| To simplify formulas, consider 
laws on No rather than N. Let us pretend for a moment that p{n) is in fact 
completely monotone, i.e.. 


r 

p{n) = / 
Jo 


^—nX 


5(A)dA 


(7.2) 


for some non-negative integrable function g{X). This holds true for the stretched- 
exponential case |19] . The probabilistic interpretation is that the law p is a mixture 
of geometric laws. Eq. (7.21 can be adopted as a definition for p(^) when Re^ > 0 
and thus provides an analytic continuation of the probability weights. Moreover we 
have, for \z\ < 1, 

5(A) 


G(z) = r 

JO 


1 — 26 


-A 


dA, 


(7.3) 


and the right-hand side provides an analytic continuation for G{z) from the unit 
disk to the slit plane C \ [1, oo). For z = e*’*'*^ with t gM. and e —>■ 0 we note 


X — t 


1 — ze 


-A 


A — t — ie 1 — e 


+ o(e)). 


(7.4) 


The function ht{X) = (A —t)/(l— ) is smooth and satisfies ht{t) = 1. Assuming 

that g is smooth—which by \7.2\ and Paley-Wiener theory directly translates into 
the decay of p(^) along lines Re^ = const—we may expect from Plemelj’s formula 
that 

lim G( e*+‘" ) = lim H ^ + ^( / (7.5) 

£\,0 £\,0 JQ X — t — !£■ X — t ) 

where the second term is the Cauchy principal value of the possibly divergent 
integral. So in particular, we expect that the limit G( e‘) = limgN^o ) exists, 

moreover 


1 

ImG(e*) = 7r5(t), p(n) = — / e "‘lmG(e‘)dt 

^ Jo 


(7.6) 


and we recover Eq. ( |4.6[ ) and some of the results from Section|^by starting from the 
relation (7.2) rather than Assumption |2 . 2| and the Lindelof integral. In view of (iTdJ 
we would like to conclude that Im G( e‘) ^ | (t) \ as in Theorem |4.4 

however doing so would require a kind of Tauberian theorem (we had worked with 
Assumption 2.3 instead). 
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In general Assumption |2. 1 1 might not be enough to guarantee a representation as 
in (7.2 1 , so following HZl one may take the approximate relation ( |7.1[ ) as a starting 
point. We introduce a truncation sequence (y„) and define 


PnW- k>y^. 


(7.7) 


Notice that YlV^iPni.^) might be slightly different from 1. By the considera¬ 
tions sketched above, for each fixed n, the associated generating function Gn{z) = 
Y^^=oPn{k)z^ extends to the slit plane C\ [l,oo), and coefficients can be expressed 
by contour integrals as we have done. The principal idea is then to show that for 
a judicious choice of (?/„) coefficients of large powers [z"^]G{z)^ behave approxi¬ 
mately like The necessary estimates have been carried out in [17] in 

a slightly different context—inversion formulas for Laplace transforms instead of 
contour integrals of generating functions—but it should be possible to adapt them 
to the present context. 

Key to the present approach is a property we might term “asymptotic complete 
monotonicity”, where the probability weights are asymptotically equivalent to the 
Laplace transform of some integrable non-negative function. This property is closely 
related to log-convexity: every completely monotone function is log-convex m 
Lemma 4.3]. The converse is, in general false [T3|, but Eq. shows that if log- 
convexity is strengthened to Assumption 2.1 then at least p{n) is asymptotically 


equivalent to a completely monotone sequence. It would be interesting to extend 
the method to more general log-convex sequences, under weaker conditions than 
Assumption |2.1[ 


7.2. Complex analysis. Another option, especially when the weights may not be 
log-convex, is to try to obtain an analytic interpolation based on the Paley-Wiener 
theorem (known as the Sampling Formula). 


Theorem 7.1. ^ If f is an entire analytic function of type < T whose restriction 
to ^ S M is then 




2 

2 



n— — oc 
oo 




sin[(^ — mr/T)T] 
^ — mr/T 

2 


n— — oo 


and 


To see how to use this, let G{z) be a probability generating function with probability 
weights p{n) = c„, n S N. Then one may define the analytic symbol for this case 
by 


^ sin[7r($ - n)\ sin[7r(^ - n)] 

^ " 7r(^ — n) 2 ” 7r(^ — n) 

n—1 ^ ' n— — oo ^ ^ 


(7.8) 


where c_„ = c„ for n 7 ^ 0 and cq = 0. Equivalently, p{f) = ^ G( e'*) e“‘^‘ dt. 
It is straightforward to check that p[f) so defined is an entire analytic function 
of type < TT. Moreover, since (c„) is square summable, the sampling theorem 
establishes that the restriction of p{f) to the real line gives an if function with 
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norm 

.00 ^00 

n——oo n—1 

In order to apply such an analytic representation to the framework outlined in 
this paper one needs to assume that Assumption 2.2(ii) holds on some appropriate 
contours in the complex plane which is needed to ensure that the Lindelof integral 
exists and has an analytic continuation to slit complex plane. The above represen¬ 
tation, in terms of entire analytic functions is probably too rigid to apply in the 
case of classes such as regularly varying hazard functions which extend the specific 
case of a sub-exponential distribution. We are currently exploring other versions of 
interpolation on restricted domains that may provide more flexibility in extending 
our results. 


Appendix A. Proofs of Lemmas 12.141 and 12.151 

Proof of Lemma \2.14\ The proof of Assumption |2.1| is straightforward and left to 
the reader. The function q{^) = is analytic in Re^ > 0 and p(^) = exp(— 
satisfies, for all fc S N, 


_ l^jfe g-|?r cos(aarg(C)) ^ |^|fc cos(a7i-/2) ^ 


(A.l) 


Since a S (0,1), we have cos ^ > 0 and Eq. (A.l) shows that is integrable 

along Re^ = 1/2 and that p(^) grows slower than any exponential exp(e|^|). This 
proves Assumption |2.2[ 

The equation q"(x*) = —l/(ncr^) can be solved explicitly. N* and N** are best 
determined with the scaling relation ( |2.24 1. They have already been determined 
in na, we omit the proof. For the insensitivity scale, we notice that 

na q [lyn) = na a (A.i) 

which goes to zero if and only if Nn S> □ 

Proof of Lemma \2.15[ The function q{x) = — logc -I- (loga;)^ is clearly smooth on 
( 1 , 00 ). Then q'{x) = /3(loga;)^“^ jx and as a: —)■ 00 , 




jSilogxY ^ 


q^^\x) 


2fi{\ogx)P 


(A.3) 


Assumption 2.1 is easily checked. For Assumption 2.2 we note q{^) = c+ (log^)^ is 
analytic in Re^ > 1. Fix 6 > 1 and write ^ = rexp(i0). As |^| —)■ 00 along Re^ = 6 
i.e. ^ = h + \y^ the argument 9 goes to ±7r/2 and we have 

Re (log 1^1 -b lef = Re (^log \y\ + ^ log(^ 

= (loglyl)'^ + 0(1), 


1 

^ —2 

yA 


+ 




(A.4) 

(A.5) 


conditions(i) and (ii) in Assumption ra are easily checked. Assumption |2.3K iii) 
follows from a computation similar to (A.3). Set yo = \Jr'^ — We have for 
i = b + iy, y > yr, uniformly in r. 


Req{^) -Req{zr) = (logy)^ - {logyof + o(l) 

> /3(log log — -b 0 ( 1 ) 
yo 


(A.6) 
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hence 


noo nc 

Jvr Ji 


’yods 


^ ^-Req{b+iyo) VO _ g-Re g(&+iyo)+0(log r) j\ 

{\ogyo)0-^ 

which proves Assumptioi] |2.3[ i) . Next let ^ G C with Re^ > 1, write ^ = rexp(i0), 
then 

iO \/3-i 


^q'iO = = /3(logr + 10 = q'(r)(l + 


and for large r and 9 G (0,7r/2), 

P-1 9 


0 


(A.8) 


(A.9) 


Im^g'(r) rlogTsinS 
and so Assumption |2.3[ ii) holds. 

We now turn to the asymptotic behavior of the sequences a;* , N* and N**. Since 
q"{xn) = it is clear that a;„ —>■ oo as n —>■ oo. The equation is 

1 /3(logx;)'5-i 




{<y 


(A.IO) 


Consequently, 

«)2 ~ Pna^{\\og{xlfy 


/3ncr^ (5 log n)^ 


log/3cr^ + (/3 — 1) loglogx* N^-i (A-11) 


1 + 


logn 


)' 


The last bracket is asymptotically equal to 1 and we get the expression for x*. 

(A.12) 


Next, we have from (2.10) 


K =< + na- 


:P{\OgXlY- 


2x1. 


The last asymptotics follows from (A.IO). 


We now turn to N**. It is asymptotically given by the solution of the equations 

(A.13) 




2na‘^ 2na^ 

t / \ A^ X. 

q (Xn) = 


'n *^71 

na"^ 


Eq. (A. 14) is equivalent to 

x^ - N„Xn + Pna-"^(logXn)'^~^ = 0. 


The relevant solution is 


i -4/3ncr2(logx„)/3-i) 

^^(loga;„)^"^(l + o(l))). 


(A.14) 

(A.15) 

(A.16) 


It follows that Nn — x„ ~ (logx„)^ We insert this in (A.13); using logx„ ^ 
log A^„, we get 

Af^ 

^ - N2(logiV„)'' - i/3'na2(logiVj2'5-2 = o(l). (A.17) 
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The relevant solution is 




na 


(logiV„)^ + V (logiV„)2/3 -h ^ 2 (log jV„)2/3-2 


r^2na^i^logN^f (A.18) 

~ 2^“^ncr^(logn + log2(T^ + /31oglogiV„)^. 

Only the term log n matters in the last bracket and the result follows. 

The last part of the lemma on insensitivity sequence is shown as in [U Section 
8.3] the proof is therefore omitted. □ 

Appendix B. Bivariate Hessian 

As explained in Step 4 of the proof outline, the Hessian at (tnjCn) has determi¬ 
nant — 1 + o(l) and is a saddle point of considered as a function of two 

real variables >0. In order to get rid of off-diagonal elements in the Hessian 
and to give all eigenvalues the same sign, we take complex ^ and reparametrize, as 
sketched in Step 5. 

Lemma B.l. Let be the unique solution o /= t — q'{^) = 0. Set 
Fn : (0, oo) X K -)■ C, Fn{t,s) = -I-is). 

Then (VF„)(t„, 0) = 0, 


HessF„(<, 0) = 


m 0 

0 q"m) 


Pit) = - 


det(Hess <l)„)(t, ^(t)) 

“ qW) ■■ 


Note that so as n —>■ oo 

detHessT’„(t„,0) =/3(t„)(?"(Cn) = - det(Hess 4>„)(f„, ^„) = l-|-o(l). 

Proof. We have 

dtFn{t,s) = (i9t$„)(t,C(t) +is) -f id^<^n)it,iit) +is)C'(t), 
dsFn{t,s) = -f is). 

At t = tn, s = 0, we have f{t) = and (VF„)(t„,0) = V<i>„(t„,^„) = 0. For the 
Hessian, we compute 

dsFnit,0) = -a|$„(t,^(t)) = q''i^it)), 


(B.l) 


(B.2) 


dtdsFnit,{)) = \{dtd^<^n)it,iit)) -f i5|$„(t,^(t)) 


5?A„(t,0) = ^$n(L^(i)) = Pit)- 


By definition of ■f(t), 


0 = nit, at)) = idtd^^n)it,f.it)) + dl^nit,f,it))i' it). 


dt 

It follows that dtdsFnit, 0) = 0, and 

12 

= idt^n)it,^it)) + idtd^<t>n)it,^it))£,'it). 


(B.3) 

(B.4) 


Pit) = ^^nit,m) = J^i^t^n)it,m) 


(B.5) 


We solve for f'(t) in Eq. (B.4), insert into Eq. (B.5), and obtain the formula for 
Pit). □ 
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